11  /  AD-AIOO  832  NATIONAL  RESEARCH  COUNCIL  OF  CANADA  OTTAM 
ff  TRAWLER  XCIN6!  A  COMPILATION  OF  WORK  DONE  A1 

II  DEC  80  J  R  STALLABRASS 

1 1  UNCLASSIFIED  DME-MD-56  NRC-19372 

ONTARIO)  --etc  F/6  A/2 
N*R.C«  (61VRA8E  DE— CTCIU 

NL 

1 

fm 

nil 

m-i 

j 

jf-l 

jra 

i 

■ 

I 

I 

I 

I 

I 

J 

I 

UL  ' 

I 

■ _ 

^  FILE  COPY 


■  ^  National  Research  Consei^hational 

■  ^  Council  Canada  \  de  recherches  C 


de  recherches  Canada 


TRAWLER  ICING 


A  COMPILATION  OF 
WORK  DONEATN.R.C. 


J.R.Stallabrass  a 
Division  of  Mechanical  Engineering 


DTIC 

ELECTEjp 
JUL  1  1981  '• 


OTTAWA 
DECEMBER  1980 


TTii*  dozumeat  has  b#en  apples 
for  puiJiic  reioci*# 
distribution  it  unlimited. 


Canada 


MECHANICAL  ENGINEERING  REPORT 
MD-56 

N.R.C.  NO.  19372  i 

$1  06  9.9  226  J 


T.R.  Ringer,  Head/Chef 

Low  Temperature  Laboratory/ 

Laboratoire  des  basses  temperatures 


E.H.  Dudgeon 
Director/Directeur 


SUMMARY 


This  report  collates  and  reviews  the  results  of  an  eleven-year 
project  on  the  icing  of  fishing  trawlers. 

Initially  the  factors  causing  the  icing  of  ships  at  sea  are  presented 
and  some  aspects  of  the  physics  of  the  icing  process  are  discussed. 

Some  simple  icing  tests  on  cylinders  of  various  diameters  are 
described.  These  were  made  in  an  icing  wind  tunnel  and  demonstrated  the 
effects  of  air  temperature  and  cylinder  size  on  the  resulting  ice  formations. 
Later,  a  number  of  possible  methods  for  reducing  the  icing  hazard  were 
tested  in  the  wind  tunnel  and  at  an  outdoor  test  rig.  Of  the  methods  tested 
the  most  effective  was  an  inflatable  rubber  de-icing  blanket. 

Throughout  the  project,  icing  report  forms  were  distributed  to 
vessels  operating  out  of  Canadian  east  coast  ports.  The  data  obtained  by  this 
means  has  helped  to  establish  the  types  of  weather  systems  responsible 
for  icing  conditions,  the  geographical  extent  of  the  occurrence  of  icing,  and 
statistics  on  the  severity  of  icing  encounters  in  the  study  area. 

Appendices  present  simplified  formulae  for  the  droplet  collection 
efficiency  of  cylinders  and  rectangular  bodies,  the  derivation  of  an  analytical 
expression  for  the  rate  of  ice  build-up  on  fishing  trawlers,  and  a  bibliography 
of  ship  icing., 


(frangais  au  verso) 


RESUME 


Dans  le  present  rapport,  nous  compilons  et  analysons  les  resultats 
d’un  projet  de  onze  ans  sur  le  givrage  des  chalutiers. 

Dans  un  premier  temps,  nous  traitons  des  facteurs  a  I’origine  du 
givrage  des  navires  en  mer  et  de  certains  aspects  de  la  physique  du  processus 
de  givrage. 


Nous  decrivons  quelques  essais  simples  portant  sur  le  givrage  de 
cylindres  de  differents  diametres.  Ces  essais  ont  ete  menes  en  soufflerie  et 
ont  permis  dt  mettre  en  evidence  les  effets  de  la  temperature  de  I’air  et  des 
dimensions  du  cylindre  sur  ia  formation  de  la  glace  resultante.  Nous  avons 
par  la  suite  effectue,  en  soufflerie  et  sur  un  greement  d’essai  place  a  I’exte- 
rieur,  des  essais  sur  un  certain  nombre  de  methodes  susceptibles  de  reduire 
le  risque  de  givrage.  Parmi  ces  methodes,  la  plus  efficace  a  etc  I’utilisation 
d’une  couverture  de  degivrage  gonflable  en  caoutchouc. 

Tout  au  long  du  projet,  nous  avons  distribue,  aux  chalutiers  cir- 
culant  au  large  des  ports  de  la  cote  et  du  pays,  des  formules  de  releve  sur  le 
givrage.  Les  donnees  ainsi  obtenues  nous  ont  aide  a  determiner  les  types  de 
systemes  meteorologiques  qui  creent  des  conditions  de  givrage,  a  evaluer 
I’etendue  geographique  ou  se  produit  le  givrage,  et  a  produire  des  statistiques 
sur  I’intensite  du  givrage  observe  dans  la  zone  etudiee. 

Les  annexes  contiennent  des  formules  simplifiees  sur  le  coefficient 
de  captation  des  gouttelettes  des  cylindres  et  des  structures  rectangulaires, 
I’eiaboration  d’une  expression  analytique  donnant  le  taux  d’accumulation 
de  glace  sur  les  chalutiers  et  une  bibliographic  des  ouvrages  portant  sur  le 
givrage  des  navires. 
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TRAWLER  ICING 

A  COMPILATION  OF  WORK  DONE  AT  N.R.C. 


1.0  INTRODUCTION 

As  a  result  of  the  loss  of  three  Canadian  fishing  vessels  ftom  icing  in  preceding  winters,  a 
research  investigation  into  ship  icing  was  proposed  in  1968  jointly  by  the  then  Department  of  Transport, 
Ship  Safety  Branch,  and  the  National  Research  Council. 

This  project  included  a  preliminary  investigation  of  the  icing  characteristics  of  cylindrical 
objects,  and  some  simple  evaluation  tests  of  a  number  of  proposed  passive  de-icing  methods.  However, 
the  main  thrust  of  this  project  has  been  the  collection  of  data  on  icing  encounters  reported  by  ships 
operating  out  of  Canadian  east  coast  ports,  with  the  objective  of  determining  the  extent  and  severity 
of  icing  conditions  in  waters  off  the  Canadian  east  coast. 

This  report  consolidates  various  reports  and  papers  written  during  the  course  of  this  project 
which  has  now  been  terminated.  A  bibliography  is  presented  in  Appendix  C. 


2.0  FACTORS  CAUSING  THE  ICING  OF  SHIPS  AT  SEA 

A  perennial  hazard  faced  by  the  crews  of  vessels  operating  in  cold  northern  waters  during  the 
winter  months  is  that  of  icing.  This  applies  particularly  to  the  Canadian  east  coast  fisherman  for,  in 
winter,  conditions  conducive  to  icing  may  exist  at  any  time  after  leaving  port,  while  steaming  to  or 
from  the  fishing  grounds  or  while  actually  fishing;  in  fact,  his  whole  trip  is  likely  to  lie  within  a  zone  of 
potential  icing. 

The  icing  of  vessels  can  result  from  a  variety  of  causes  as  indicated  by  a  number  of  authors 
(Refs.  1  to  6)  as  follows: 

(a)  Supercooled  fog  (commonly  referred  to  as  arctic  frost  smoke  or  black  frost). 

(b)  Freezing  rain  or  freezing  drizzle. 

(c)  Falling  snow,  in  particuleir  wet  snow. 

(d)  Freezing  sea  spray  (vessel  and/or  wind  generated). 

Of  these,  the  freezing  of  sea  spray  on  a  vessel  is  the  most  dangerous  cause  of  icing.  The  other 
three  causes,  all  of  which  result  in  fresh  water  ice,  are  generally  not  considered  seriou.s.  However, 
several  writers,  particularly  Borisenkov  (Ref.  3)  and  Shekhtman  (Ref.  8),  indicate  that  the  most  serious 
type  of  icing  occurs  when  one  or  more  of  these  three  lesser  cau.ses  occurs  in  combination  with  freezing 
sea  spray.  On  the  other  hand,  it  has  been  suggested  by  Tabata  (Ref.  2)  that  snow  and  fog  should  not  he 
emphasized  in  connection  with  the  icing  of  ships. 

There  is  general  agreement  regarding  the  specific  meteorological  factors  that  result  in  the 
freezing  of  sea  spray  on  a  ship  (Refs.  3,  4,  7,  8).  The  air  temperature  must  he  below  the  freezing 
temperature  of  the  sea  water  (about  -1.8°C,  but  dependent  on  .salinity),  the  sea  temperature  lower 
than  about  +6°C  (this  limit  is  less  well  defined  and  depends  on  other  factors,  particularly  air  tempera¬ 
ture),  and  the  wind  speed  greater  than  about  20  knots  (Beaufort  Force  5),  although  at  lower  wind 
speeds  some  light  icing  is  likely  to  occur. 

In  addition  to  these  factors,  which  apply  also  to  ocean  structures,  the  degree  of  icing  s('verity 
experienced  by  a  vessel  is  also  dependent  on  the  speed  and  heading  of  the  vessel  relative  to  tlie  wind  as 
well  as  on  the  size  and  lines  of  the  vessel  (including  freeboard,  amount  of  rigging  and  otht'r  top  hampi'r, 
etc.). 


-2- 


Combined  with  the  relative  wind  speed  and  direction  factors,  the  wave  height  and  fre()uency 
also  have  a  bearing  on  the  severity  of  icing,  since  these  together  determine  the  amount  and  distribution 
of  water  and  spray  breaking  over  the  vc'ssel.  Steep  waves  of  short  period  result  in  more  spray  than  hmg- 
period  waves  or  swell. 

It  is  apparent  from  all  this  that,  with  so  many  variables,  many  of  which  are  not  easily  meas¬ 
ured  or  characterized,  the  prediction  of  the  rate  of  icing  on  a  vessel  is  no  simple  task. 

3.0  THE  ICING  PROCESS 

The  formation  of  ice  on  a  body  results  when  water  impinges  on  the  body  in  a  sub-freezing 
environment  and  the  rate  of  loss  of  heat  from  the  body  is  such  that  some  or  all  of  the  water  is  frozen 
before  it  can  run  off  the  body.  Clearly,  the  colder  the  water,  the  greater  the  chance  that  this  will  occur. 
If  the  impinging  water  is  supercooled  (i.e.  its  temperature  is  below  its  equilibrium  freezing  point)  then 
it  is  a  certainty  that  ice  will  be  formed  providing  there  is  no  heat  input  to  the  process,  such  as  provided 
by  a  heated  body  or  by  the  kinetic  heating  of  high  speed  impact  between  the  water  and  the  liody.  These 
cases  need  not  concern  us  in  a  discussion  of  the  icing  process  (although  a  liody  might  be  heated  as  an 
icing  prevention  method). 

The  sea  water  may  be  several  degrees  above  its  equilibrium  freezing  iioint.  but  due  to  wave 
and  wind  action  and  to  the  (lassage  of  the  vessel,  water  is  detached  and  some  of  it  strikes  the  vessel. 

The  decks  may  be  washed  by  quantities  of  "green”  water  which  has  cooled  little  and  may  or  may  not 
create  much  ice.  .Vs  height  increases,  the  quantity  of  water  becomes  less  and  its  degree  of  suiniivision 
becomes  greater,  with  the  result  that  much  of  this  water  will  be  cooled  substantially,  and  th('  smaller 
the  drop  size  the  more  nearly  will  its  temperature  approach  that  of  the  ambient  air.  These  drops  are 
likely  to  remain  in  a  supercooled  liquid  state*  during  their  flight  time  from  parent  sea  to  striking  the 
vessel's  structure.  In  general  then,  the  greater  the  height  above  the  sea.  th(‘  lowa-r  the  water  concentra¬ 
tion  will  be.  but  a  greater  proportion  of  that  water  will  freeze  and  form  ice. 

3.1  Drop  Deposition 

As  the  size  of  the  water  drops  decreases,  aerodynamic  drag  forces  increase  relative  to  tlie 
inertia  forces  acting  on  the  drops  and  there  is  a  tendency  for  the  trajectories  of  the  drops  to  more 
nearly  conform  to  the  streamlines  of  the  air  as  it  is  deflected  around  an  object.  The  larger  the  object, 
the  greater  the  deflection  for  a  given  dro[)  size  and  v(>locity  and  the  lower  the  proportion  of  droiis  tliat 
impinge  on  the  object.  Thus,  the  actual  impingement  or  catch  rate  is  less  than  the  potential  catch  rate 
assuming  straight  line  motion;  the  ratio  of  the  actual  to  potential  catch  is  known  as  the  catch  efficiency 
or  collection  efficiency.  Figure  1  shows  the  streamlines  of  air  around  a  cylinder  and  the  manner  in 
which  the  droplet  trajectories  are  deflected  from  a  straight  path.  Only  those  droplets  lying  between 
the  tangent  trajectories  impinge  on  the  cylinder,  thus  the  Collection  Efficiency,  E,„  ,  is  defined  by  the 

S 
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As  the  drops  coo!  below  their  equilibrium  freezing'  temperature  they  do  not  immediately  tree/e  but  exist  in  a  metastahle  super¬ 
cooled  liquid  state.  Pure  water  can  support  a  maximum  superco<dinp  ot'about  4(r('  tor  small  droplets  ol  a  lew  microns  tli.irncfcf 
and  about  33'^(*  for  volumes  of  a  lew  cubic  centimetres  I  Ref.  9).  However,  in  nature  ft>reien  particles  in  the  water  act  as  free/inc 
nuclei  causing  free/ine  to  occur  at  temperatures  much  warmer  than  those  for  homoeenei>us  nuclcation.  nevertheless  sipnilicant 
supercooling  still  occurs,  the  most  likely  frce/in^  temperature  being  in  inverse  propv>rtion  to  the  logarithm  ol  ilroplct  si/c  I  lie 
effect  of  salts  in  solution  on  the  other  hand  is  to  lower  the  nucleation  leinperalure  (  Ref.  IP).  )  lie  etieci  ol  lliese  tw <'  oppc'smc 
phenomena  on  the  supercooling  and  nucleation  ol  sea  water  drops  is  unknown  but  is  obvunislv  ot  significance  to  the  siud>  ot  sea 
spray  icing.  It  has  been  suggested  that  below  abt>ut  ^18''('  the  spray  is  frozen  and  strikes  the  ship  as  small.  ilr>  crystals  that  do  mit 
adhere  ( Ref.  1  1 ),  however,  this  is  ct>ntradicled  by  Shekbiman  I  Kel.  8J  who  cites  instances  of  icing  occurring  at  air  temperatures 
down  to 
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The  equations  of  motion  of  particles  in  a  stream  of  air  appear  to  have  been  first  developed 
by  Albrecht  (lief.  12)  and  later  applied  by  Taylor  (Ref.  13)  to  droplets  causing  icing  of  aircraft.  Values 
of  E,„  are  not  easily  calculated;  however,  Langmuir  (Ref.  14)  presented  droplet  trajectory  results  in 

graphical  form  for  simple  objects  such  as  spheres,  cylinders,  and  flat  ribbons  in  terms  of  two  non- 
dimensional  parameters.  The  use  of  graphical  results  of  this  sort  is  not  always  convenient,  particularly 
in  numerical  computations,  so  that  in  Appendix  A  (taken  from  Ref.  15)  empirical  e.xpressions  are 
presented  that  permit  the  direct  computation  of  the  collection  efficiency  of  cylinders  and  flat  surfaces 
over  a  range  of  sizes  and  conditions  applicable  to  the  icing  of  ships  and  stationary  structures,  and  witli 
sufficient  accuracy  in  view  of  present  uncertainties  with  regard  to  defining  the  spray  parameters. 


3.2  Heat  Balance 

When  the  supercooled  drops  strike  the  cold  obstacle,  it  acts  as  a  nucleating  agent  to  initiate 
freezing,  so  that,  depending  on  the  initial  .supercooling  of  the  impinging  water  and  the  various  heat 
losses  and  gains  at  the  icing  surface  (Refs.  16,17)  some  or  all  of  the  impinging  water  freezes  in  its 
immediate  area  of  impingement.  The  proportion  that  does  actually  freeze  on  impact  is  called  the 
freezing  fraction,  n,  while  the  remaining  fraction  (1-n)  flows  away  from  its  impingement  zone  under 
the  influence  of  gravity  or  of  aerodynamic  forces;  some  of  this  run-off  will  flow  over  areas  of  greater 
heat  lo.ss  where  further  freezing  will  occur  while  some  may  actually  flow  or  be  blown  completely  off 
the  surface  subject  to  impingement. 

The  freezing  fraction  and  hence  the  rate  of  icing  on  the  obstacle  is  given  by  tht'  solution  ot 
the  heat  balance  equation,  which,  if  only  the  primary  heat  transfer  i>rocesses  acting  at  the  icing  surface 
are  taken  into  consideration  (i.e.  ignoring  heat  transfer  due  to  radiation  and  kinetic  and  viscous  heating 
which  are  negligible  in  this  application),  may  be  expressed  by: 

Ri  +  +  q.  +  Re  +  q.,  =  0  '1' 

where  q,  =  latent  heat  of  freezing  a  certain  fraction,  n,  of  the  impinging  water. 

q^^  =  heating  (or  cooling)  of  the  impinging  water  to  the  equilibrium  surface 
temperature, 

q^.  =  heat  transfer  by  convection  to  the  surrounding  air, 
q^,  =  heat  transfer  by  evaporation  to  the  surrounding  air, 

q_i  =  heat  conducted  to  or  from  the  icing  surface  through  the  underlying  structure. 

The  fraction  of  the  impinging  water  that  freezes  decreases  both  as  the  temperature  increases 
and  as  the  water  concentration  increases.  The  excess  water  will  run  off  or  be  blown  off  the  surface,  and 
the  mass  of  ice  forming  on  the  obstacle  in  time,  t,  is 

M,  =  nM,^  =  nE,,,  VAwt  (2) 

where  n  =  freezing  fraction 

=  mass  of  water  impinging  on  obstacle  in  time,  t 
E|„  =  collection  efficiency 

V  =  velocity  of  air  relative  to  obstacle 

A  =  projected  frontal  area  of  obstacle 

w  =  mass  of  water  droplets  in  unit  volume  of  air  (liquid  water  content). 


As  the  ice  forms  on  the  obstacle,  the  shape  presented  to  the  flow  of  droplets  is  modified, 
and  the  projected  frontal  area  may  increase  in  size.  Thus  in  Equation  (2),  the  quantities  n,  ,  and  A 

may  all  change  with  time,  so  that  the  instantaneous  icing  rate,  Rj  =  dMj/dt,  also  becomes  a  function 
of  time. 

If,  in  Equation  (2),  A  is  considered  a  constant,  equal  to  the  frontal  area  of  the  ice-free 
obstacle,  then  it  is  possible  for  E^,  and  even  the  Icing  Efficiency,  N  =  nE,„  ,  to  assume  values  greater 

than  unity  if  the  fronted  area  of  the  iced  obstacle  exceeds  the  ice-free  area  A. 

The  above  discussion  has  looked  at  the  gross  or  overall  icing  rate  on  an  object,  and  where,  in 
Equation  (2),  the  factors  n  and  E,„  are  the  effective  or  total  freezing  fractions  and  collection  efficiency 

for  the  object  as  a  whole.  No  information  is  forthcoming  from  such  considerations  on  the  local  growth 
rate  of  ice  on  different  parts  of  the  object,  and  hence  on  the  overall  shape  of  the  ice  deposit,  except  in  a 
rather  empirical  way  (Ref.  18). 

Determination  of  local  growth  rates  requires  a  knowledge  of  the  factors  relating  to  the  local 
convective  heat  transfer  and  to  the  local  collection  efficiency.  Such  factors  are  complex  and  are  far 
from  understood,  except  perhaps  for  the  simple  case  of  a  cylinder  for  which  an  attempt  (Ref.  19)  has 
been  made  to  numerically  model  the  local  icing  rate  and  so  predict  ice  shapes  on  a  cylinder. 

4.0  ICING  TESTS  ON  CYLINDERS 

As  an  introduction  to  the  National  Research  Council’s  programme  of  ship  icing  research,  a 
series  of  tests  was  made  in  an  icing  wind  tunnel  (1.37  m  square  test  section)  to  study  the  ice  formation 
and  rate  of  growth  on  cylinders  of  various  diameters.  These  tests  are  described  fully  in  Reference  17. 
and  their  salient  details  are  presented  here. 

Five  cylinders,  3.8  cm  (IV2  in.),  7.6  cm  (3  in.),  15.2  cm  (6  in.),  30.5  cm  (12  in.),  and  45.7  cm 
(18  in.)  in  diameter  were  used,  and  one-hour  tests  were  run  on  each  cylinder  in  lioth  a  horizontal  and 
vertical  orientation.  A  wind  speed  of  43  knots  (Beaufort  Force  9)  was  used  throughout  and  each  test 
was  repeated  at  two  temperatures  (about  -  14°C  and  about  -  7°C).  The  water  sprays  were  set  to  give 
a  water  concentration  in  the  air  of  3.2  g/m^  with  an  effective  drop  diameter  of  about  0.2  mm.  In  the 
light  of  the  paucity  of  quantitative  data  on  sea  spray  conditions,  these  values  were  chosen  rather  in¬ 
tuitively  within  the  capabilities  of  the  test  facility,  but  were  considered  to  result  in  fairly  representative 
icing  conditions.  One  area  where  conditions  where  clearly  not  representative  was  that  fresh  water  was 
used  rather  than  sea  water. 

Typical  ice  deposits  are  shown  in  Figures  2  and  3,  and  representative  shapes  are  sketched  in 
Figure  4.  The  effect  of  greater  blow-off  and  run-off  of  water  at  higher  temperatures  is  evident  in  a 
smaller  ice  accretion  than  at  lower  temperature,  and  the  effect  of  gravity  results  in  an  asymmetrical 
formation  and  icicles  on  the  horizontal  cylinders,  and  a  variation  of  thickness  with  height  on  the 
vertical  cylinders. 

In  analyzing  the  results,  the  concept  of  Icing  Efficiency  was  used  (see  Section  3.2  above),  but 
in  the  form  of  a  one-hour  icing  efficiency,  defined  as: 


where 


mass  of  ice  formed  in  one  hour. 


The  ice  deposits  wee  weighed  at  the  end  of  each  1-hour  run  and  the  1-hour  icing  efficiency 
determined.  In  spite  of  the  difference  in  the  shape  of  the  ice  accretions  on  horizontal  and  vertical 


-5- 


cylinders  (as  a  result  of  gravity  affecting  the  run-off),  no  significant  difference  in  icing  efficiency  was 
apparent.  It  was  therefore  possible  to  combine  the  results  from  both  orientations  as  in  Figures  5  and  6, 
which  show  the  effect  of  temperature  and  of  cylinder  diameter  on  the  1-hour  icing  efficiency. 


No  great  significance  should  be  read  into  the  actual  rates  of  icing  of  the  test  cylinders,  since 
their  relationship  to  rates  of  icing  at  sea  is  unknown.  However,  the  icing  efficiencies  derived  from  these 
tests,  and  presented  in  Figures  5  and  6,  are  considered  to  give  a  useful  guide  to  the  fraction  of  water  that 
actually  freezes  on  rails,  masts,  lines,  etc.  over  an  extended  period  of  time.  As  noted  in  Reference  20, 
“there  is  virtually  no  knowledge  as  to  the  rate  at  which  spray  strikes  the  exposed  surfaces  of  a  ship  or  of 
the  percentage  that  freezes”.  These  remarks  still  apply  some  20  years  later;  however,  these  results 
provide  some  indication  of  the  percentage  of  spray  that  may  freeze. 


The  manner  in  which  the  icing  efficiency  increases  with  decreasing  cylinder  size,  even  to 
values  greatly  in  exeess  of  100  percent  for  small  diameter  cylinders  (Fig.  6),  demonstrates  very  force¬ 
fully  the  absolute  necessity  for  keeping  the  number  of  rails,  stays,  etc.  to  a  minimum,  and  that,  where 
it  is  not  possible  to  dispense  with  them  completely,  one  larger  diameter  strut  should  be  made  to  do  the 
job  of  several  smaller  ones.  It  was  found  that  a  12-fold  increase  in  the  basic  cylinder  diameter  from 
3.8  cm  to  45.7  cm  resulted  in  only  a  3-fold  increase  in  the  mass  of  ice  collected. 

It  was  noticeable  that  the  ice  thickness  was  not  significantly  affected  by  the  cylinder  diameter, 
but  only  by  the  temperature.  On  the  other  hand,  the  width  of  the  ice  accretion  (dimension  A,  Fig.  4) 
was  affected  by  cylinder  diameter  as  well  as  temperature,  and  it  is  this  dependence  on  cylinder  diameter 
that  is  reflected  in  the  icing  efficiency.  Figure  7  demonstrates  the  effects  of  temperature  on  ice  thick¬ 
ness,  and  of  cylinder  diameter  on  the  ice  width  expressed  as  the  ratio  of  width  to  cylinder  diameter. 

The  fact  that  the  ice  thickness  appears  to  be  independent  of  cylinder  diameter  (a  relationship 
that  does  not  apply  when  all  the  water  freezes  in  its  area  of  impingement)  leads  to  the  speculation  that, 
for  sea  spray  ice,  the  rate  of  ice  growth  might  give  an  approximate  measure  of  the  spray  concentration 
providing  the  air  temperature  and  velocity  are  known.  Tests  at  other  water  concentrations  and  speeds 
are  required  to  ascertain  the  reliability  of  such  a  measure. 

The  main  findings  and  conclusions  from  these  tests  are: 

1.  The  dependence  of  the  icing  efficiency  on  the  cylinder  diameter  and  on  the  air  temperature 
was  demonstrated;  the  icing  efficiency  being  shown  to  increase  with  both  decreasing  dia¬ 
meter  and  decreasing  temperature. 

2.  It  was  shown,  for  the  conditions  of  the  tests,  that  the  ice  thickness  is  a  function  of  the 
temperature,  but  is  independent  of  the  cylinder  diameter. 

3.  The  width  of  the  ice  accretion  is  a  function  of  both  temperature  and  cylinder  diameter,  such 
that  for  small  diameter  cylinders  the  ice  width  may  grow  to  many  times  the  cylinder  diameter. 

4.  The  results  clearly  indicate  the  necessity  for  reducing  rigging,  spars,  etc.  to  a  minimum  as  one 
very  necessary  measure  for  diminishing  the  hazard  of  icing  at  sea. 


5.0  METHODS  OF  REDUCING  THE  ICING  HAZARD 

One  method  of  reducing  the  icing  hazard  was  demonstrated  in  the  tests  just  described,  that  is 
that  the  number  of  small  scale  collectors  such  as  rigging  and  railings,  and  top  hamper  generally,  should 
be  reduced  to  a  minimum  and  where  complete  elimination  is  not  feasible  one  larger  scale  component 
should  be  made  to  serve  the  function  of  several  small  scale  ones.  This  has  been  recognized  for  a  con¬ 
siderable  time  (Ref.  11),  but  that  it  has  not  always  been  fully  taken  to  heart  would  seem  to  be  evident 
from  some  of  the  discussion  in  the  final  report  of  a  British  inquiry  into  trawler  safety  (Ref.  21). 


Unless  a  vessel  can  completely  avoid  icint;  conditions,  some  means  of  comhattintj  the  ice  are 
necessary.  Methods  may  be  subdivided  into  two  categories,  de-icing  and  anti-icing.  C^ertain  de-icing 
methods  require  the  active  intervention  of  crew  members  on  deck:  these  include  such  methods  as  the 
use  of  steam  or  water  lioses  as  well  as  the  traditional  use  of  mallets  and  a.\es.  Conditions  on  deck  at 
these  times  are  particularly  hazardous,  and  l)eeome  completely  impossible  wlien  icing  conditions  are 
the  most  severe  and  de-icing  is  the  most  necessary.  Also  such  methods  are  not  amenable  to  the  de-icing 
of  the  upper  parts  of  superstructure  and  rigging  where  the  added  weight  of  ice  contributes  most  to  tin* 
loss  of  stability.  Preferred  methods  therefore  should  preclude  the  use  of  crew  on  deck  and  so  must  be 
of  an  automatic  or  semi-automatic  nature.  Recent  proposals  or  developments  include  the  use  of  plastic 
foam  as  a  passive  de-icer  ( Refs.  22,  23),  inflatable  rubber  de-icer  blankets  ( Ref.  21 ),  and  an  electrical 
impulse  method  (Ref.  25).  The  cost  factor  seems  to  militate  against  the  extensive  use  of  electrical 
heating  (Ref.  26)  for  de-icing. 

Tests  at  NRC  (described  in  Ref.  27  and  summarized  below)  have  shown  the  remarkable 
effectiveness  of  the  inflatable  rubber  de-icer  on  both  cylindrical  and  flat  surfaces.  Little  is  known  about 
the  Russian  electro-impulse  system,  but  if  the  claims  made  for  it  are  valid,  it  might  well  he  effective  for 
the  de-icing  of  large  flat  metal  surfaces. 

Anti-icing  methods,  in  which  ice  is  prevenUxI  from  forming,  may  be  either  chemii'al  or 
thermal.  Chemical  methods  imply  basically  the  use  of  freezing  jtoint  depressant  chemicals  ait|)lied  at 
the  surface  to  be  protected,  although  they  may  be  extended  to  include  materials  to  which  icc  has 
negligible  adhesive  strength.  The  major  problem  in  the  use  of  freezing  point  de|iressants  (i.c.  glycols 
and  alcohols)  is  that  of  uniformly  distributing  the  fluid  over  the  surface.  Materials  promoting  low  ice 
adhesion  are  theoretically  of  two  types,  fluid  and  solid,  the  fluid  tyjie  being  typically  a  silicone  oil 
which  prevents  ice  from  adhering  to  the  underlying  solid  substrate;  such  a  method  although  resulting 
in  extremely  low  adhesion  (Refs.  28,  29,  30)  is  perhaps  only  practical  for  moderately  small  mechanical 
devices  on  which  an  oil  film  can  be  maintained  .so  as  to  prevent  the  mechanism  from  being  jammi'd  liy 
ice,  A  solid  material  to  which  ice  will  not  adhere  is  a  dream  not  yet  realized,  and  probably  unri>aiizal>lc. 
Certain  fluoro-carbon  polymers  do  display  a  reduced  ice  adhesion  of  aliout  15  kN/m‘  (20  jisi)  in  shc'ar. 
but  this  is  still  far  in  excess  of  the  approximate  1  kN/m^  necessiU'y  to  promote  ready  ice  release. 

Certain  flexible  silicone  coatings  hav('  shown  low  values  (about  3.5  kN/m’ )  on  initial  ice  release,  but 
this  deteriorates  on  subsequent  releases  and,  in  addition,  problems  of  ultraviolet  degratlation  and  low 
abrasion  resistance  currently  preclude  these  materials  from  (iractical  consideration. 

Thermal  anti-icing  methods  are  perhaps  rather  more  attractive  than  chemical  ones.  Large 
.scale  electrical  heating  is  perhaps  prohibitively  expensive;  however,  for  small  items  it  may  [irovide 
certain  conveniences,  a  pcU-ticular  ca.se  in  point  being  the  use  of  an  electrically  conducting  transpiumit 
layer  deposited  on  glass  and  used  for  heating  bridge  windows  (Ref.  31,  Fig,  11 ).  It  .seems  likely  that 
waste  heat  from  the  vessels'  main  engines  could  be  employed  to  good  imriiose,  either  to  heat  com¬ 
ponents  such  as  the  mast,  etc.  directly,  or  perhaps  indirectly  by  such  means  as  a  convective  thermo¬ 
syphon  as  described  by  Lock  (Ref.  32),  or  a  two-phase  thermosyphon  (Ref.  33),  the  action  of  which 
is  demonstrated  in  Figure  8.  The  two-phase  thermosyphon  has  the  advantage  of  higher  heat  transfer 
rates  and  of  remote  siting  of  the  heat  exchanger  and  boiler  section  as  long  as  the  vapour  transiiort  and 
liquid  return  are  unobstructed.  The  working  fluid  could  be  ammonia,  but  R21  would  be  more  suited 
to  the  temperatures  involved  since  the  pressures  in  the  system  would  then  be  clo.se  to  atmospheric  and 
pressure  vessel  requirements  could  be  less  stringent.  Tests  at  NRC  have  demonstrated  the  anti-icing 
effectiveness  of  such  a  system  when  applied  to  simulated  masts  and  handrails. 

5.1  Tests  Made  at  NRC 

Evaluation  tests  (described  in  detail  in  Ref.  27)  were  made  on  a  number  of  mi'thods  proposed 
for  reducing  the  hazard  of  superstructure  icing  on  shijjs  at  sea.  The  tests  were  mad(‘  under  simulated 
icing  conditions,  either  in  the  1.37  m  icing  wind  tunnel,  or  at  an  outdoor  test  site  using  natural  low 
temperatures  in  the  wintertime.  This  site  is  shown  in  Figure  9.  Tmi  [lanels  1.2  m  X  0.9  m  were  attached 
to  a  vertical  frame  of  approximately  2.5  m  X  4.6m.  Various  methods  of  de-ii’ing  were  mounted  on  or 
applied  to  the  panels.  Two  cables  (a  parallel  filament  rope  and  a  stei'l  cable)  were  strung  diagonally 


forward  and  down  from  the  two  upper  corners  of  the  frame.  To  augment  the  prevailing  wind  a  fan  was 
used  and  on  its  front  face  spray  nozzles  were  mounted  to  produce  the  icing  spray.  In  both  the  wind 
tunnel  and  outdoor  facility,  fresh  (tap)  water  was  used  to  produce  the  icing  conditions. 

The  methods  tested  were  all  de-icing  methods,  that  is  to  say  they  were  all  concerned  with  the 
removal  of  the  ice  after  it  had  formed  on  various  structures,  rather  than  with  the  prevention  of  its 
formation.  The  specific  methods  employed  were: 

1.  Pneumatic  de-icer  as  applied  to  a  mast 

2.  Plastic  foam  as  applied  to  a  mast 

3.  Plastic  foam  with  rubber  surface  as  applied  to  a  mast 

4.  Pneumatic  de-icer  as  applied  to  a  bulkhead,  such  as  a  bridge  front 

5.  Plastic  foam  as  applied  to  a  bulkhead 

6.  Plastic  foam  with  rubber  surface  as  applied  to  a  bulkhead 

7.  Paint  and  varnish  on  a  wooden  bulkhead 

8.  The  distribution  of  a  freezing  point  depressant  (ethylene  glycol)  over  a  bulkhead 

9.  Parallel  filament  rope  for  use  as  stays. 

5.1.1  Pneumatic  De-icer 

A  pneumatic  de-icer  is  essentially  a  tube  or  series  of  tubes  built  into  a  rubber  mat,  such  that 
when  the  tubes  are  inflated  with  air  they  ex[iand  and  the  mat  surface  heaves  and  stretches,  and  in  so 
doing  breaks  down  the  adhesion  of  the  ice  to  the  mat  surface.  Such  pneumatic  de-icer  mats  have  long 
been  used  for  de-icing  the  leading  edges  of  aircraft  wings,  and  recently  their  application  has  been  vs- 
tended  to  marine  use.  In  this  application,  their  use  has  been  proposed  for  the  de-icing  of  masts  and 
stays,  bridge  fronts,  radar  scanners,  and  liferaft  stowages  (see,  for  instance.  Ref,  34  or  35). 

These  tests  constitute  an  independent  evaluation  of  their  effectiveness  in  de-icing  masts  and 
flat  surfaces.  The  test  specimen  was  a  short  length  of  mast  de-icer,  122  cm  long  and  91  cm  wide, 
suitable  for  attachment  to  a  30.5  cm  diameter  mast.  Details  are  given  in  Figure  10.  The  same  de-icer 
was  used  both  for  mast  de-icing  tests  in  the  wind  tunnel  and,  when  mounted  on  a  flat  surface,  for 
bridge  front  de-icing  tests  at  the  outdoor  test  stand. 

Si.\  tests  were  made  in  the  wind  tunnel  on  the  pneumatic  de-icer  installed  on  a  30.5  cm  dia¬ 
meter  cylinder,  137  cm  long,  to  represent  a  short  length  of  mast.  The  de-icer  was  not  bondcKi  to  the 
cylinder,  liut  was  laced  onto  it  in  the  manner  shown  in  Figure  11.  The  cylinder  was  mounli'd  in  the 
tunnel  in  such  a  way  that  it  could  be  oscillated  about  its  normal  i)osition  through  ±180°;  in  this  way 
ice  could  be  built  completely  airound  the  mast,  a  situation  that  might  occur  should  a  vessel  change 
course  a  number  of  times  during  an  icing  encounter.  Thus  it  was  possil)le  to  see  whether  there  was  any 
change  in  the  de-icer’s  performance  between  complete  and  incomplete  encapsulation  of  the  mast  in  ice. 
Quite  low  temperatures  (-9  to  -23°C)  and  a  high  water  concentration  (3.2  g/m’ )  were  u.sed  for  thest' 
tests  in  order  that  the  de-icer  be  subjected  to  a  severe  test. 

Only  two  of  the  wind  tunnel  mast  tests  will  be  described  here.  In  the  first,  at  23°(’,  the 
mast  was  exposed  to  the  icing  condition  for  nearly  two  hours  during  which  lime  it  was  rotated  *  180 
to  produce  an  ice  coating  uniformly  tliick  around  the  whole  circumference.  The  radial  tliickm-ss  of 
the  ice  at  mid-height  was  about  13  cm  on  the  peaks  and  10  cm  in  the  valleys  of  the  rough  ice  formation 
(Fig.  12).  A  further  five  minutes  were  allowed  to  elap.se  after  the  sprays  were  cut  off  to  allow  all  water 
on  the  surface  to  freeze  and  some  cooling  of  the  ice  to  occur.  One  inflation  of  the  de-icer  with  an  air 
pressure  of  100  kPa  (15  psi)  broke  this  thickness  of  fully  encapsulating  ice  cleanly  and  decisively,  as 
Figure  12  shows. 

The  second  test  to  be  described  was  at  a  high  temperature,  10°C.',  and  for  a  duration  of 
only  45  minutes.  Because  of  the  higher  temperature,  the  water  flowed  on  the  ice  surface  to  a  greater 
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extent  than  on  the  earlier  test.  The  mast  was  not  rotated  for  this  test  and  so  a  rather  flat-faced  accretion 
resulted  (Fig.  13),  being  about  4  cm  thick  at  the  centre  front,  but  increasing  to  9  cm  thickness  at  45°  on 
either  side,  and  tapering  off  to  nothing  at  about  ±100°.  An  air  pressure  of  70  kPa  on  the  de-icer  shed  this 
ice  cleanly  with  a  single  inflation. 

Eight  tests  were  made  at  the  outdoor  test  site  on  the  pneumatic  de-icer  configured  as  a  bulk¬ 
head  de-icer.  The  de-icer  was  mounted  flat  on  one  of  the  1.2m  X  0.9m  panels  by  clamping  to  the  panel 
the  bars  that  passed  through  the  attachment  ferrules  down  each  side  of  the  de-icer.  The  de-icer  was  not 
bonded  to  the  panel  so  as  to  permit  the  required  deformation  when  the  inflation  air  pressure  was  applied. 

Because  of  the  variability  of  wind  direction,  intensity,  and  gustiness,  even  when  augmented  by 
the  fan,  it  was  not  possible  to  assess  an  effective  wind  velocity  for  each  run;  it  was,  however,  generally 
in  the  order  of  10  m/s  (20  kt).  For  this  reason,  and  because  throughout  these  tests  alterations  were  made 
in  the  water  spray  system  to  improve  the  spray  distribution  over  the  test  panels,  the  spray  concentration 
was  highly  variable  from  run  to  run.  Since  these  tests  were  very  qualitative  in  nature,  it  was  felt  that 
these  deficiencies  were  not  of  great  consequence  as  long  as  the  ice  formation  produced  appeared  reason¬ 
ably  characteristic  of  shipboard  icing. 

Again  only  representative  tests  will  be  described.  In  one,  the  de-icer  was  mounted  on  the 
upper  centre  panel.  About  10  cm  of  ice  was  allowed  to  form  at  a  temperature  of  -  15°C,  and  a  long 
overnight  cold  soak  at  -  19°C  was  allowed  before  de-icing  was  attempted.  Thus  an  abnormally  severe 
test  was  given  the  de-icer,  requiring  more  than  one  inflation.  The  first  inflation  cracked  the  ice  on  the 
de-icer,  but  it  remained  anchored  around  the  edges  and  did  not  shed.  A  second  inflation  caused  a  small 
amount  of  ice  to  break  away  (Fig.  14b,  c)  but  not  until  the  third  inflation  was  the  bulk  of  the  ice  dis¬ 
lodged  (Fig.  14d).  Since  the  ice  was  anchored  at  the  edges,  and  particularly  at  the  top  where  the 
de-icer  tended  to  sag  away  from  the  steel  panel  behind,  there  seems  little  doubt,  had  the  de-icer  been 
much  larger  and  had  it  been  subjected  to  the  buffeting  experienced  at  sea,  that  the  ice  would  have 
shed  on  the  first  inflation. 

In  other  tests,  two  pneumatic  de-icers  were  used,  one  either  side  of  the  upper  centre  panel 
onto  which  was  bonded  rubber-faced  plastic  foam.  It  was  thought  that  the  simultaneous  inflation  of 
the  two  de-icers  might  also  remove  the  ice  from  the  panel  between.  Three  tests  were  made  with  this 
arrangement  at  de-icing  temperatures  of  -12°C,  -8°C  and  -3°C.  At  -  12°C,  ice  4  cm  thick  was  shed 
cleanly  from  one  de-icer  while  ice  on  the  other  de-icer  remained  anchored  along  the  top  and  sides, 
although  being  detached  from  the  active  surface  of  the  de-icer  (Fig.  15).  At  -8°C  ice  shed  cleanly  from 
both  de-icers,  but  only  at  -3°C  when  the  ice  was  rather  wet  and  rotten  was  the  desired  result  of  removing 
the  ice  from  the  intervening  rubber-covered  panel  achieved. 

5.1.2  Plastic  Foam 

As  a  result  of  Japanese  sea-going  tests,  mats  of  plastic  foam  applied  to  the  hull,  bridge  front, 
breakwaters,  deck,  etc.  are  reported  to  be  “effective  in  the  prevention  of  icing,  and  furthermore  make 
ice  removal  extremely  easy”  (Ref.  22). 

It  was  known,  before  Reference  22  was  available,  that  foamed  plastic  was  being  experimenUnl 
with  in  Japan,  but  it  was  not  clear  exactly  how  it  functioned  as  an  ice  release  mechanism,  nor  the 
precise  materials  u.sod.  However,  some  foamed  polyethylene  sheet  was  obtained  for  ex|)eriinental  use. 
This  material  was  10  mm  thick  and  had  a  specific  gravity  of  0.05. 

This  plastic  foam  .sheeting  was  tested  both  in  the  wind  tunnel  on  a  30  cm  diameter  steel 
cylinder  representing  a  mast,  and  on  flat  steel  plates  on  the  outdoor  test  stand.  In  both  instances  the 
foam  was  bonded  onto  the  steel  substrate  using  a  rubber-based  adhesive.  T(‘sts  were  made  both  on  the 
bare  plastic  foam  and  with  a  0.8  mm  thick  sheet  of  neoprene  rubber  cemented  to  the  outer  surface. 

Two  test  runs  were  made  in  the  wind  tunnel  on  the  bare  plastic  foam  at  a  temperature  of 
about  14°C.  Other  conditions  were  as  for  the  pneumatic  de-icer.  The  cylinder  was  not  rolatfxl  and 
the  ice  grew  to  a  thickness  of  3  cm  on  the  front  centreline  in  20  minutes  of  icing.  As  soon  as  the  tunnel 
was  shut  down  after  Run  1,  the  ice  was  assailed  with  a  wooden  club.  The  ice  was  rather  mu.shy,  and 


-9  - 

although  it  tended  to  be  compressed  somewhat  by  clubbing,  it  was  removed  moderately  easily.  After 
Run  2,  20  minutes  were  allowed  for  the  ice  to  harden  before  removal  was  attempted.  This  time  the  ice 
broke  off  in  chunks  without  much  effort;  however,  no  large  area  of  ice  was  removed  in  one  piece, 
making  removal  rather  tedious. 

Next,  the  polyethylene  foam  was  covered  with  the  sheet  of  neoprene  rubber,  which  gave  a 
smoother  and  more  durable  surface  to  the  foam.  Two  more  icing  runs  were  made.  Thirty  minutes  were 
allowed  after  Run  3  for  the  ice  to  harden.  The  effort  required  to  remove  the  ice  was  judged  to  l)e  about 
the  same  as  for  Run  2.  Only  five  minutes  were  allowed  after  Run  4,  which  was  made  at  a  temperature 
of  -8°C,  and  the  ice  removal  was  judged  much  easier  than  after  Run  1. 

In  the  tests  at  the  outdoor  test  site  it  was  found  that  ice  could  only  be  removed  with  some 
difficulty  from  polyethylene  foam  bonded  to  a  flat  surface,  although  with  less  difficulty  than  from 
plain  painted  steel  panels.  The  compressible  foam  beneath  the  ice  layer  permitted  the  ice  to  be  cracked 
and  broken  with  a  club  or  hammer  more  readily  than  on  an  inflexible  steel  panel;  however,  the  broken 
pieces  of  ice  generally  remained  adhered  to  the  rather  rough,  accommodating  surface  of  the  foam,  and 
it  was  usually  necessary,  after  hammering,  to  pull  the  broken  pieces  off  the  foam  laboriously  by  hand. 

In  addition,  the  shattering  of  the  ice  frequently  cut  the  foam. 

Ice  removal  from  the  rubber-surfaced  foam  was  somewhat  less  difficult  and  showwi  a  decided 
temperature  dependence.  After  a  run  at  about  -  19°C  it  was  noted  that  the  5  cm  thick  hard  porous  ice 
was  moderately  difficult  to  remove  using  a  wooden  club,  while  ice  on  the  steel  panels  was  practically 
impossible  to  break  and  remove  by  this  method.  Ice  was  finally  removed  by  a  sledge  hammer.  Kven 
then,  only  small  amounts  of  ice  could  be  broken  away  with  each  blow.  The  ice  was  judged  .slightly 
easier  to  remove  at  about  -  16°C,  although  thin  ice  (about  1  cm)  just  splintered,  leaving  pieces  still 
adhered  to  the  rubber. 

At  -8°C  ice  could  be  removed  without  too  much  trouble;  after  initially  breaking  it  with  the 
club,  the  ice  was  lifted  off  easily  by  hand.  At  -  3°C  the  ice  was  very  soft  and  wet  and  was  easily  rt'moved 
from  all  panels. 


5. 1 .3  Paint  and  Varnish 

The  basic  finish  on  the  flat  steel  panels  at  the  outdoor  test  site  was  a  black  ant'-rust  paint. 
During  the  course  of  the  tests  one  of  the  panels  was  repainted  with  a  grey  d(>ck  f»ainl.  and  woo<len 
panel  was  finished  half  in  the  grey  deck  paint  and  half  in  varnish. 

There  was  no  discernible  difference  in  the  removal  of  ice  from  any  of  the  steel  panels,  whether 
painted  with  the  black  anti-rust  paint  or  the  grey  deck  paint  they  were  all  equally  difficult.  Remov:il 
became  easier  as  the  temperature  increased,  and  little  difficulty  was  experienced  at  3^('  m  removing  in' 
from  any  of  the  panels. 

The  ice  was  not  quite  as  difficult  to  remove  from  the  wooden  iianel,  but  this  was  thought  to 
be  due  to  its  lower  flexural  rigidity  compared  with  that  of  the  steel  panels.  If  was  judgeil  that  the  ice 
released  slightly  more  easily  from  the  grey  paintiHl  half  of  th<' {)anel  than  from  the  varnished  lialf. 
although  the  difference  was  so  slight  as  to  be  of  no  real  significance. 

5.1.4  Freezing  Point  Depressant 

Kthylene  glycol  was  distributed  over  one  area  of  the  vertical  steel  panel  wall  of  the  ovitdoor 
test  stand.  The  glycol  was  fed  under  gravity  from  a  tank  |)laced  80  cm  above  a  horizontal  distribution 
pipe  from  which  it  was  discharged  through  a  series  of  1  mm  diameter  holes  spaced  at  T.O  cm.  The 
discharge  holes  were  so  oriented  that  the  glycol  flowed  onto  the  panel  .surface,  and  as  it  flowed  ilow  n 
the  panel  surface  the  individual  streams  from  each  hole  l«*nded  to  spread  and  join,  particularly  when 
the  surface  was  wet. 


-  10- 


In  theory  the  freezing  point  depressant  should  mix  completely  with  the  impinging  water  and 
form  a  mixture  of  which  the  freezing  point  is  lower  than  the  ambient  temperature,  thus  resulting  in 
complete  freedom  from  icing.  In  practice  it  is  not  possible  to  ensure  ideal  distribution,  quantity,  or 
mixing  of  the  depressant.  Thus  ice  will  form,  but  this  ice  is  usually  infested  with  the  depressant  and  so 
is  mushy  and  has  poor  adhesion  to  the  surface  on  which  it  forms;  removal  is  therefore  easy. 

This  was  found  to  be  the  case.  In  one  run  at  '8°C,  about  a  quarter  of  the  panel  was  free  of 
ice  at  the  end  of  the  test,  which  deposited  about  2.5  cm  of  ice  over  the  test  area.  Elsewhere  on  the  panel 
the  ice  was  ju.st  standing  clear  of  the  panel  and  only  required  a  touch  to  collapse  it.  In  another  test  at 
10°(’  glycol  was  fed  at  a  rate  of  1.1  litre/hr  to  six  panels,  i.e.  an  average  rate  of  1  litre/hr  per  6  m’ .  .At 
the  end  of  the  run,  the  glycol-infested  ice  was  broken  away  from  the  panels  with  great  ease. 

5.1.5  Parallel  Filament  Rope  and  Steel  Cable 

A  parallel  filament  rope  was  compared  with  a  steel  cable.  7'his  plastic  rope  consists  of  a 
load-bearing  core  of  highly  compressed  parallel  filaments  of  polyethylene  terephthalate  (a  linear  poly¬ 
ester)  covered  by  a  black  polyethylene  sheath.  To  quote  the  manufacturer,  “A  twistless  rope  has  been 
evolved  which  is  light,  strong,  flexible,  maintenance-free,  and  most  important  —  highly  resistant  to 
stretch”.  Most  significant  from  the  icing  point  of  vh'w  was  the  smooth  polyethylene  sheath  and  its 
flexibility,  which  suggested  the  possibility  of  low  adhesion  between  ice  and  rope. 

A  6  m  length  of  2-ton  rope  (11  mm  0.1).)  was  .strung  from  the  top  right-hand  corner  of  the 
mounting  frame  for  the  wall  ()anels  to  the  front  UTt  corner  of  the  decking  (directions  as  viewed  from 
the  front),  while  for  compariscm  a  6  m  length  of  12.7  mm,  6  X  19  galvanized  steel  cable  was  strung 
from  the  tof)  left  corner  to  the  front  right  corner  (Fig.  15). 

Because  of  the  difference  in  torsional  ngkiity  l)etween  the  parallel  filament  rope  and  the 
steel  cahlt'.  some  difference  in  the  manner  of  icing  was  observetl.  The  rope,  because  of  its  lowiT  tor¬ 
sional  stiffiK'ss.  lend<“d  to  twist  under  the  w<‘ight  of  ice  forming  on  its  windward  side,  and  in  so  tloing 
more  of  its  circumferi'nce  would  become  coated  will)  ice  than  was  the*  case  with  the  sti-el  cable  whose 
stiffness  [irecluded  significant  twisting  in  llie  short  length  used.  Tliis  twisting  effi'ct  was  more  niiuked 
at  lower  tenqieratures  or  lower  water  concentrations  (Fig.  11)  under  wliich  conditions  mori'  of  the 
im))inging  water  froze  on  the  windward  side  than  at  higher  temiii'ratures  and  higher  water  concentra¬ 
tions.  when  most  of  the  ice  that  formi'd  on  the  calilc's  was  in  thi'  form  of  pendant  icick's  that  excTted 
little  torsional  moment  (Fig.  15). 

In  extreme  cases  complete'  circumferc-ntial  enca|)sulation  of  the  [ilaslic  rope-  occurred  as  a 
result  of  the  twisting.  In  one  test  ciicatisulation  occurred  tcjward  the*  lower  end  of  the'  ro|)e  .so  that 
when  it  was  givc-n  a  good  shake  by  hand  only  the  ice  on  the  uppc'r  part  of  the-  rope  was  di.slodged. 

Since  the*  rope  and  the  steel  calile  contacted  at  their  mid  points,  shaking  the  ro(M‘  also  .shook  the  cable, 
with  the*  result  that  most  of  the  ice  on  the  cable  was  also  .shed.  The  ici'  that  shi'd  had  left  the  roju' 
com|>letely  clean  and  fn-e  of  ice  remnants,  indicating  true*  adhesional  rek'a.se,  while  small  ice  remnants 
still  adhered  to  the  steel  cable,  indicating  that  the  ice  had  failed  cohesively  rather  than  adhesivc'ly  at  the 
ice/metal  interface.  In  this  particular  test,  ice  on  the  lower  j>arl  of  the  steel  calde  was  uj)  to  15  cm  thick, 
with  icicles  15  or  more  cm  long. 

It  was  found  on  this  and  other  occasions  that,  if  ice  at  the  lower  end  of  the  parallel  filament 
rope  was  removed  by  ra[)|)ing  the  icetl  ro{)e  with  a  club,  encapsulating  ice  further  ui>  th<'  ro)M'  would 
slide  down  and  could  in  turn  be  bri  ken  off.  On  nt)  occasion  was  the  statii-  weight  of  the  ice  sufficient 
to  ove-'-ome  the  adhesive  force  lietween  the  ice  and  the  rope,  ami  so  promote  a  self-protective  action. 


5.1,6  Discussion  and  Conclusions 

The  evaluation  of  these  various  m<'thods  was  very  qualitative  and  the  test  conditions  may  not 
have  resembled  actual  sea  icing  conditions  with  any  degree  of  exactness.  However,  com|)arative  evalua¬ 
tions  are  thought  to  be  rea.sonably  valid.  The  use  of  fresh  water  instead  of  .salt  water  is  thought  to  have 
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made  the  tests  somewhat  more  severe  by  virtue  of  the  greater  strength  of  fresh  water  ice.  Other 
mitigating  circumstances  at  sea  are  the  vibration  and  buffeting  experienced  in  rough  weather  —  condi¬ 
tions  that  could  not  be  simulated  on  the  simple  outdoor  test  facility. 


The  pneumatic  de-icer  proved  most  effective,  both  as  a  mast  de-icer  and  on  flat  surfaces. 

There  is  little  doubt  that  the  occasional  failure  to  shed  ice  on  the  first  inflation  of  the  de-icer  would  not 
have  occurred  under  the  buffeting  of  actual  sea-going  conditions  and  with  a  larger  active  area  than  used 
in  these  tests.  Disadvantages  of  this  otherwise  most  effective  method  of  de-icing  are;  1)  cost,  2)  the 
likelihood  of  damage  to  the  de-icers  by  cutting  and  abrasion  if  used  in  or  adjacent  to  work  areas,  3 )  the 
danger  to  hands  working  on  deck  from  ice  falling  from  above,  and  4)  the  need  for  removal  of  the  shed 
ice  from  the  decks  if  reduction  of  the  total  ice  weight  is  to  be  achieved. 

On  the  other  hand,  pneumatic  de-icing  appears  to  be  an  ideal  method  of  de-icing  the  upper 
parts  of  masts,  etc.  to  reduce  the  large  heeling  moments  caused  by  ice  on  these  parts. 

The  use  of  ethylene  glycol  as  a  freezing  point  depressant  proved  effective  in  reducing  the 
strength  and  adhesion  of  ice  on  flat  surfaces.  It  is  probable  that,  under  the  buffeting  of  sea-going 
conditions,  a  fair  measure  of  anti-icing  performance  would  be  achieved  liy  this  method  with  the  flow 
rates  used  in  these  tests.  This  method  of  icing  prevention  may  prove  useful  for  small  components  such 
as  radomes  and  inflatable  life  raft  containers.  Its  disadvantages  when  used  in  work  areas  are:  1 )  a  slip¬ 
pery  residue  could  he  left  on  the  decks  after  an  icing  encounter  if  not  completely  flushed  away  by  water, 
and  2)  a  danger  exists  of  contamination  of  the  fish  catch. 

The  various  coatings  and  finishes  tested  were  not  considered  very  successful.  The  most  suc¬ 
cessful  was  the  rubber-covered  plastic  foam.  This  was  becau.se  the  deformable  foam  sulistrate  p(“rmitt(‘d 
breaking  of  the  ice  with  moderate  ease,  while  the  ice  exhibited  lower  adhesion  to  lh('  smooth  rubber 
surface  than  to  the  pel)bled  surface  of  the  bare  polyethylene  foam.  Th(‘  rubber  also  provided  a  durable 
surface  to  the  otherwise  easily  damaged  foam.  No  case  occurretl  of  ice  falling  off  the  foam  or  rul)ber- 
covered  foam  under  its  own  weight,  as  suggested  by  the  Japanese  work;  this  may  be  due  partly  to  the 
small  panel  size,  partly  l)ecause  of  th('  lack  of  vibration  and  buffeting,  and  also  perhaps  bc'causi*  )nost 
of  the  tests  were  made  at  temperatures  well  below  the  fn'ezing  point  with  fresh  wat(“r  instead  of  sea 
water. 


In  addition,  although  on  the  Japanese  vessels  of  steel  construction  tin-  mats  were  cementf'd 
on,  in  the  case  of  the  vessels  of  wooden  construction  the  mats  were  affixed  with  steel  straps  at  tin*  top. 
bottom,  and  sides  of  each  mat,  and  the  resulting  diaphragm  action  of  these  mats  under  the  action  of  the 
sea  undoubtedly  was  the  rea.son  that  “there  was  little  icing  on  the  mats,  and  ice  was  sprung  off  every 
time  there  was  any  kind  of  shock  or  shaking”. 

A  disadvantage  of  this  method  is  the  problem,  also  encountered  in  the  Japanese  tests,  of 
effectively  bonding  or  attaching  the  foam  and  rubber  blankets  to  flat  surfaces.  However,  litth'  difficulty 
is  seen  in  lacing  blankets  on  masts,  etc.,  much  in  the  manner  used  for  the  pneumatic  de-icer.  This  might 
|)rove  u.seful  in  work  areas  where  the  use  of  pneumatic  de-icers  may  be  prohibited  becau.se  of  the  damage 
hazards. 


If  the  various  methods  are  rated  in  order  of  de-icing  effectivene.ss  or  ease  of  ice  removal,  the 
following  list  results: 

1.  Pneumatic  de-icer 

2.  Freezing  point  depressant 

3.  Rubber-surfaces  pla.stie  foam  on  steel  panel 

4.  F’aint  or  varnish  on  wooden  panel 
Bare  polyethylene  foam  on  steel  panel 

6.  Paint  on  steel  panel. 


-12- 


In  spite  of  the  obviously  lower  adhesive  strength  of  ice  to  the  polyethylene  sheath  of  the 
parallel  filament  rope  than  to  the  zinc  galvanizing  of  the  steel  cable,  little  difference  was  observed  in 
the  ease  of  removal  of  ice,  even  though  it  released  cleanly  from  the  rope,  while  cohesive  failure  of  the 
ice  resulted  in  residual  ice  fragments  remaining  on  the  steel  cable.  During  these  tests,  encapsulation 
occurred  on  the  plastic  rope  only,  because  of  its  lower  torsional  rigidity;  however,  in  stays  that  are 
more  vertical,  twisting  would  be  less  likely,  but  ice  encapsulation  might  be  more  likely  to  occur  as  a 
result  of  the  flow  of  water  around  the  cable.  Under  these  circumstances  difficulty  would  certainly  be 
encountered  in  removing  ice  from  steel  cables;  however,  the  tendency  for  ice  to  slide  down  the  rope 
once  ice  at  a  lower  level  had  been  removed  would  permit  the  removal  of  ice  from  otherwise  out-of-reach 
locations.  A  major  disadvantage  of  the  parallel  filament  rope  is  the  ease  with  which  it  can  be  cut,  thus 
limiting  the  locations  at  which  it  could  be  employed  aboard  ship. 


6.0  LOCAL  CONDITIONS  OFF  EASTERN  CANADA 

The  sea  area  of  interest  in  this  study  lies  approximately  between  the  latitudes  42° N  and  60° N 
and  between  the  45°W  meridian  and  coast  of  North  America.  In  this  area  the  combination  of  strong 
winds  and  low  air  temperatures,  two  of  the  factors  required  for  icing,  are  provided  by  the  passage  of 
intense  cyclones.  Archibald  (Ref.  36)  states  that,  for  a  five-year  period  (1963-67)  for  January,  February 
and  March,  the  predominant  path  of  the  most  intense  storms  lies  right  along  the  Nova  Scotia  coast  and 
across  the  southern  part  of  Newfoundland,  with  a  somewhat  less  frequent  track  lying  parallel  and  south¬ 
eastward  of  this  predominant  track.  This  is  illustrated  in  Figure  16  which  is  reproduced  from  Archibald’s 
paper.  The  dots  with  concentric  circles  around  them  indicate  the  position  of  maximum  intensity  of  the 
centre  of  the  most  intense  storms,  i.e.  those  with  a  pressure  at  the  centre  of  between  956  and  972  mb 
and  with  a  very  steep  pressure  gradient.  Less  intense  storms,  with  pressures  of  between  972  and  980  mb, 
are  represented  by  solid  dots.  It  is  clear  that  many  of  these  storm  centres  will  produce  in  their  wake 
strong  northerly  or  northwesterly  winds  carrying  cold  Arctic  air  to  the  sea  area  of  interest.  The  result  is, 
therefore,  that  these  regions  are  particularly  susceptible  to  the  meteorological  conditions  most  conducive 
to  icing. 


The  oceanographic  conditions  of  interest  are  the  sea-surface  temperature,  the  sea  state  and,  to 
a  lesser  extent,  the  sea-surface  salinity. 

Monthly  mean  sea  surface  temperatures  are  presented  in  Reference  37.  Those  for  Febniary, 
which  do  not  differ  very  significantly  from  those  for  January  and  March,  are  shown  in  Figure  17.  These 
show  cold  water,  with  temperatures  down  to  about  -1°C  adjacent  to  land  or  to  the  sea  ice  cover,  slowly 
increasing  in  temperature  with  distance  from  land  such  that  the  8°C  isotherm  lies  about  150  nautical 
miles  off  the  southern  tip  of  Nova  Scotia  (Cape  Sable),  some  250  nautical  miles  off  Cape  Race,  New¬ 
foundland,  and  considerably  more  removed  to  the  east  of  Newfoundland  and  Labrador. 

The  roughness  of  the  sea  depends  not  only  on  the  wind  strength,  but  also  on  the  fetch  and 
the  duration  of  the  wind  (Ref.  38).  This  is  illustrated  in  Table  I.  Some  sheltering  of  the  fishing  grounds 
off  Nova  Scotia  and  south  of  Newfoundland  from  northerly  winds  is  afforded  by  the  proximity  of  the 
land;  a  large  part  of  the  winter  fishing  is  done  within  100  nautical  miles  of  land,  resulting  in  significant 
reduction  in  potential  wave  height  for  winds  of  Force  6  and  above.  However,  a  large  part  of  the 
St.  Pierre  Bank  (45°  —  47°N,  55°  —  57°W)  and  the  Scatari  Bank  (46°N,  59°W)  is  not  protected  from 
the  north-westerly  winds  which  sweep  across  the  Gulf  of  St.  Lawrence  and  through  the  Cabot  Strait. 
Similarly  the  fetch  may  be  considered  unlimited  over  most  of  the  Grand  Banks  and  for  sea  areas  east  of 
Newfoundland. 

The  freezing  temperature  of  sea  water  is  a  function  of  its  salinity  and  may  be  reprc'sented  by 
the  expression : 


t  =  -0.002  -  0.0524S  0.00006S^ 


where  S 
and  t 


salinity  in  parts  per  thousand  (%<■), 
equilibrium  freezing  temperature  (°C). 
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This  expression  is  within  0.002°C  of  the  tabulated  values  of  the  freezing  (melting)  tempera¬ 
ture  of  sea  water  as  quoted  by  List  (Ref.  39)  for  salinities  from  0  to  40%i>. 

The  mean  isohalines  (Ref.  37)  for  the  three  winter  months  January  to  March  indicate  a  con¬ 
siderable  variation  in  the  area  of  interest,  from  about  30%o  (parts  per  thousand)  in  the  Cabot  Strait  to 
about  35%o  at  the  southern  edge  of  the  Grand  Banks.  However,  this  variation  in  salinity  only  results  in 
a  variation  in  the  freezing  temperature  of  about  a  quarter  of  a  degree,  i.e.  from  about  -  1.65°C  to  -  1.9°C, 
as  Figure  18  demonstrates,  hardly  significant  enough  to  have  any  effect  on  the  onset  or  the  severity  of 
icing. 

These  considerations  lead  to  the  conclusion  that  conditions  conducive  to  frequent  and  severe 
icing  occur  in  the  area  under  consideration  during  the  winter  months. 

To  investigate  these  conditions  further,  to  compile  statistics  on  the  frequency  and  severity 
of  icing  encounters,  and  to  determine  actual  geographical  extents  of  icing  condictions,  data  have  been 
collected  by  means  of  voluntary  reports  using  a  standard  report  form.  The  next  section  will  discuss 
this  programme. 


7.0  ICING  REPORTING  PROGRAMME 

Icing  report  forms  were  distributed  jointly  by  the  Ship  Safety  Branch  of  the  Canadian  Coast 
Guard  and  the  Division  of  Mechanical  Engineering  of  the  National  Research  Council  of  Canada  to  all 
vessels,  but  primarily  fishing  vessels,  operating  out  of  Canadian  east  coast  ports.  Forms  were  distributed 
to  some  350  fishing  and  other  vessels  over  20m  in  length.  Of  these,  about  240  are  30m  or  longer,  and 
it  is  estimated  about  70%  may  be  at  sea  at  any  one  time.  Reporting  was  voluntary,  and  as  a  result  rather 
variable,  both  as  regards  compliance  and  the  quality  of  the  response;  however,  in  general,  both  the 
quantity  and  the  quality  of  the  responses  improved  with  time.  Reports  were  received  from  only  about 
15%  to  20%  of  the  vessels  over  30m  length,  while  on  any  one  day  no  more  than  about  10%  of  the  vessels 
actually  at  sea  made  out  reports,  even  when  icing  conditions  were  severe. 

The  report  forms  themselves  underwent  several  modifications  during  the  survey  period  in  an 
attempt  to  make  them  simpler  and  more  acceptable  to  the  user,  wiiile  at  the  same  time  providing  as 
much  relevant  data  as  possible.  The  current  version  of  the  form  is  shown  in  Figure  19.  Earlier  versions 
were  presented  in  Refs.  15  and  40. 

Ships’  masters  were  urged  to  make  returns  even  when  no  icing  was  encountereii  on  a  tri)! 

(i.e.  a  null  report)  so  that  a  better  appreciation  of  the  frequency  of  icing  encounters  could  be  gained. 

In  general  only  a  handful  of  masters  were  assiduous  in  this,  so  that  frequencies  of  occurrences  are  based 
only  on  a  proportion  of  the  vessels  responding,  the  number  of  which  varied  from  year  to  year,  as  ilid  the 
area  in  which  the  vessels  operated. 

Table  II  shows  the  number  of  vessels  responding  each  year  as  well  as  the  number  of  reports 
submitted. 

7.1  Severity  Statistics 

Based  on  the  reports  from  all  vessels  for  .s«>ven  winters  and  on  tin-  reported  si-vc-rity  of  i<’ing 
(i.e.  the  subjective  a.s.sessment  given  by  the  ship’s  officer  making  the  report),  statistics  on  tlu-  daily 
expectation  of  various  degrees  of  assessed  severity  for  each  of  the  wint<-r  months  have  been  compilt-<l. 
The  results  are  presented  in  Table  III.  For  each  of  the  months,  December  to  March,  and  for  the  com¬ 
plete  4-month  period,  1  December  to  31  March,  the  Table  gives  the  number  of  days  that  icing  of 
severity  light  or  greater,  moderate  or  greater,  or  heavy  was  reported.  The  yearly  average  number  of 
days  and  standard  deviation  from  this  average  is  shown.  The  great  variability  from  year  to  year  for  tin- 
months  of  December  and  March  is  evident  by  the  high  standard  deviations  (approximately  90'’;  of  tin- 
averages).  Much  greater  yearly  consistency  is  evident  for  the  months  of  January  and  February,  Based  on 
on  these  7  years  (1970-71,  1971  72,  1973-74,  1975-76,  1976-77,  1977-78,  and  1978-79)  tin-  daily 
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expectation  of  experiencing  icing  of  the  three  severity  categories  is  given;  for  example,  in  December 
more  than  trace  icing  may  be  expected  on  about  15%  of  the  days  (i.e.  about  5  days)  while  in  February 
this  expectation  is  some  50C4  (14  days).  Similarly,  for  the  4  month  season,  the  daily  expectation  is  36'> 
(43  days)  for  light  icing  or  greater,  while  heavy  icing  may  be  expected  on  11  days  (9'^  ). 

It  should  be  noted  that  the  data  of  Table  [II  relate  to  icing  reported  from  throughout  the 
study  area.  Because  vessels  are  not  necessarily  in  every  part  of  the  study  area  at  all  times,  these  statistics 
may  in  some  cases  rather  underestimate  the  frequency  of  occurrence. 

Frequencies  in  specific  areas  within  the  total  study  areas  may  be  expected  to  be  lower  than 
the  overall  frequencies  given  in  Table  III.  Two  specific  areas  have  Ijeen  selected  for  similar  analysis, 
namely  Divisions  3L  and  3P  of  the  area  designated  by  the  International  Commi.ssion  for  the  Northwest 
.Xtlaiitic  Fisheries  (ICNAF),  and  shown  in  Figure  22.  Taltle  IV'  jiresents  the  results  for  the  K'.N’.VF  Div¬ 
ision  3L  and  suggests  frequencies  of  about  one  third  of  those  for  the  total  study  area,  with  frecpiem  ies 
of  about  IS'^r  (4  to  5  days  per  month)  for  icing  of  at  least  light  severity  during  the  months  of  January 
and  February.  In  February  the  frequency  drops  to  5.6''!  for  moderate  or  greater  severity  with  large 
yearly  variability,  the  likely  cause  of  this  being  the  encroachment  of  .sea  ice  into  this  area  to  a  greater 
extent  some  years  than  others.  The  results  for  the  ICNAF  Division  31’  (the  .Sub-divisions  3I’n  and  3l’s 
are  combined)  are  given  in  Table  V.  In  this  sea  area  the  fre(juencies  of  occurrence  are  seen  to  be  gener¬ 
ally  slightly  over  one  half  of  those  for  the  total  study  area,  except  for  the  month  of  D('cember  when 
frequencies  of  less  than  2%  are  ex|)erienced.  The  higher  frequencies  for  December  seen  in  Table  111  for 
the  total  area  result  largely  from  icing  reports  from  more  northerly  s(‘a  areas. 

The  foregoing  statistics  relate  to  the  subjective  icing  severity  assessments  as  madi'  by  tlie 
ships’  officers,  and  there  is  little  reason  to  question  tlucse  assessments  in  vitnv  of  th(>  great  expi'rit’iici' 
these  men  have  of  icing  conditions  each  winter  and  of  the  ice  handling  qualities  of  their  rc'spective 
vessels.  However,  these  statistics  relate  more  to  the  rate  of  icing  experienced  ratlier  than  to  the  tottil 
ice  accumulated  on  the  vessels.  .Since  the  amount  of  accumuhited  ice  is  tlie  significant  factor  its  hir  its 
the  stability  and  handling  qualities  of  the  vessel  are  concernetl,  and  is  of  special  significance  to  reguhtting 
bodies  concerned  with  the  safthy  of  ships  at  sea,  statistics  Itased  on  ice  accumulation  hitvt*  biH'u  compiieti. 

Since  most  nations  with  significant  dee|)  sea  fishing  fleets  subscribe  to  the  recommcniialions 
on  the  stability  of  fishing  vessels  of  the  Inter-Oovernmental  Maritime  Consultative  Organization  tlMCOl. 
it  has  been  foutid  convenient  to  express  the  degr(*e  of  icing  in  terms  of  the  ’MC'(4  icing  allowance  crite¬ 
rion.* 


Bast'd  on  seven  years  of  data  using  returns  from  only  those  vessels  for  which  a  relativeix 
complete  history’  of  voyages  throughout  the  winter  sea.son  was  available.  Table  \'l  has  been  drawn  up 
anti  uptlates  that  presentetl  in  Reference  42.  The  tiegree  of  icing  is  ba.setl  on  the  thicknesses  of  ice  on 
various  parts  of  the  vessel  as  reportetl  in  the  tiuestionnaire  (Fig.  19).  and  converletl  rather  arbitrarily  to 
multiples  or  sub-multi|)les  of  the  IMCD  criterion  using  the  points  system  illustratetl  in  Figure  20.  The 
api>lication  of  this  points  system  is  most  easily  tlenninstrated  by  the  use  of  an  example.  In  this  i-xample 
it  will  be  assumed  that  the  ftillowing  ice  thicknes.si's  have  been  reported: 

a)  Diameter  of  ice  on  forward  rails  T.b-lOO  mm  (3-,b  in.) 

b)  Diameter  of  ice  on  other  rails  T.b-lOO  mm  (3-5  in.) 


1  111-  Kcci'iniiK’nddt HHi  «mi  Iiii.k  t  Sl.ihililv  «*l  !  isluni:  W’vscK  <  Kcl .  4  1 1  .tf  the  Inlcr^  lOvciniiK'ni.i!  M.inlinu’  (  <'n''iilld!  i\  lOiiMin/.i- 
fit  HI  I  !M<  O)  lav  s  tliiwn  flic  fiillms  me  fiiinimum  icc  .illov^.incc  m  cakuLitmc  a  vcwcl'v  nIaIhUIv  . 

•i)  kc  in^  oil  cxpt'scil  wc.ifhcr  ticiks  .iml  vMni'w.jys. 

h )  >  kf  m  '  ft >r  flic  pri'ict  IctI  l.itcr.il  .mM  »»f  c.k  h  suK  of  the  vessel  .ibi»v e  the  \v .iler  pl.iiic 


(  Mferion  .i)  .ih<*ve  implies  .m  .ivcr.iee  it  c  l  hick  ness  .»!  .ibt>uf  in  in  •*!  I  'v  in  w  Inie  t  t‘irespontfs  f>'  .tii  .iv  ci.h’c  iIik  kiu  ss 
of  rttuehlv  mm  or  H  m 
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c)  Diameter  of  ice  on  forward  stay 

d)  ’riiiekness  on  main  deek 

e)  'I'hiekness  oti  boat  lieek 

f)  TIriekness  on  wheelhouse  font 
tJ)  'I’hiekness  on  Imlwarks 


>  100  mm  (>  5  in.) 
50-75  mm  (2-3  in.) 
25-50  mm  (1-2  in.) 
50-75  mm  (2-3  in.) 
25-50  mm  ( 1-2  in.) 


Thus,  rising'  Figure  20,  (a)  is  rated  I  [mints,  (b)  8  points,  (e)  8  points,  (d)  8  [toints,  (e)  4  (mints, 
(1)8  ()oints  and  (R)  4  (mints.  'I’lie  [loints  total  is  tberc'fore  4-1,  wliieli  falls  in  the  ranfte  3()-48,  thus  m- 
ttieating  a  rating  ecpiivalent  to  1'  >  times  the  IMC’O  ieing  allowanee. 

Fable  VI  shows  the  percentage  of  voyages  on  whieh  ieing  of  the  indicated  degree  was  en¬ 
countered  by  the  assemblage  of  vessels.  If  ieing  was  eneounti'red  more  llian  onei’  on  a  voyage,  only  the 
more  severe  encounter  was  used  in  these  .statieties.  'Fhi'  [x'riod  covered  is  a  ty()ieal  ieing  season,  from 
roughly  the  middle  of  December  to  the  etui  of  .March  (a[)()ro.\imately  100  liays).  During  this  time  an 
individual  vesst'l  makes  some  7  or  8  Iritis  having  a  tytiical  duration  of  about  10  to  12  days. 

The  totals  for  the  seven  years  for  whieh  statistics  could  be  obtained  show  there  is  a  40% 
(irobability  of  icing  of  some  degree  being  encountered  on  a  voyage  to  the  fishing  grounds  during  the 
icing  sea.son.  It  is  also  shown  that  on  8%  of  winter  voyages  tlu>  dc'gree  of  icing  encountered  will  e.xcecd 
th('  INK'O  it'ing  allowance,  while  2.5'^  of  voyages  will  ex()erienc<'  icing  of  greater  than  twice  the  IMC'O 
standard. 


7.2  (leographical  Distribution  of  Icing 

'File  locations  of  icing  encountt'rs  given  by  .sonu>  150  .shi()  icing  retiorts  are  shown  in  Figuia'  21. 
'Fhey  indicate  that  icing  may  be  encountered  anywh(>re  within  the  coniiiu'ntal  shelf  /.one  of  Kastern 
Canada  as  well  as  the  adjacent  <l(‘e|)er  wat(>rs  of  the  Labrador  Sea.  'Fhe  very  large  sea  t('m|)eratur(>  gra¬ 
dients  existing  off  the  Scotian  Slnlf  and  the  .Ntnvfoundland  Grand  Banks  as  a  n>sult  of  the  warm  Gulf 
Stream  (irovide  a  very  ilefinite  limit  to  the  southern  extent  of  icing  (Hef.  5). 

For  more  (h'tailed  study,  the  s(>a  ar(*as  in  which  icing  occurrences  have  been  re()orted  are 
conveniently  defined  in  terms  of  the  Divisions  designated  by  the  Inti'rnational  Commission  for  the 
Northwest  .Atlantic  Fisheries  (ICN.XF).  Figure  22  is  a  ma[)  showing  tlu'.st'  Divisions.  In  ri'|)orting  the 
areas  of  ii-ing  occurrence's  it  is  found  convenient  to  combine  Divisions  31’n  and  31’s  into  a  single  Dj'.  i- 
sion  3P,  and  to  redefine  the  boundary  between  Divisions  IVn  and  IV’s  by  moving  it  south  to  the  (larallel 
■l  l’lO'N,  suc'h  that  the  bulk  of  the  continental  shelf  within  4V  is  in  the'  redc'fined  Division  IVn. 

The  ri'tmrts  suggest  that  there  are  (ireferred  fishing  areas,  which  nc'verthc'less  tend  to  vary 
somewhat  from  year  to  year  and  even  from  month  to  month,  [lartly  as  a  result  of  the  availability  of 
fish  stocks  ;ind  also  as  a  result  of  the-  encroac  hment  of  sc'a  ice,  both  .southwccrcl  along  the  coasts  of 
Labrador  and  Newfoundland  and  within  the-  Gulf  of  St.  Lawrc'ncc'.  'Fhus  it  is  not  (lossible  to  inter(irc't 
the  total  number  of  icing  rc'tiorts  from  various  arc-as  as  indicating  the'  relative  frec|uency  and  or  severitv- 
of  icing  in  the-  res()ectivc'  areas  at  any  given  time*,  ilowc’ver,  by  taking  retiorts  of  icing  in  inclivuluai 
ICN.-\F  Divisions  for  each  of  the  wintc'r  months  ovc'r  a  numbe-r  of  yc'ars  it  becomes  inissible  to  deter¬ 
mine'  whether  Ihc're  is  a  likelihood  of  ext>c'rieneing  icing  in  a  givc-n  Division  and  month.  Fight  years  of 
data  liave  lic'c'ii  analy/eci  by  Division  and  month  to  dc'terminc'  whc'ther  any  casc-s  of  icing  of  severitv 
grc'ater  than  slight  (i.e.  half  IMCO  criterion  or  grc'atc'r)  have  bc'c'ii  rc'tiortc'd.  Fhe  results  arc-  (iresc'iited 
in  Table'  VII  and  illustratc'd  as  (lercc'iitagc's  for  each  month  in  Figure's  23  to  28. 

Since'  a  reejuirc'd  condition  for  se-rious  icing  according  to  various  source's  ( Hi'fs.  3.  4,  7.  8)  is 
a  combination  of  air  tc'mtie'raturi'  below  the'  fre'czing  (>oiiu  of  sc'a  watc'r  (about  1  .S'C)  and  se'a-surface 
te'mtie'raturc'  below  about  -t-ri  ’(',  it  is  to  In'  ex(n'i't('ei  that  such  a  combination  could  also  In'  usi'd  to  de'- 
fini'  arc'as  of  e'X(n'cte'd  icing.  Finis  in  Figure's  29  to  ,31  isothe'rms  of  40°F  ( 4.4  ’(')  nu'an  monthly  st'a- 
surface'  ti'nitii'raturi'  as  givc'ii  in  Ki'fc're'nce'  37.  and  iso()leths  of  10';  (n'rce'ntage'  fre'e|ui'ncy  of  air 
ti'mtn'rature's  of  O’C  or  Ic'ss  as  give'll  in  Iti'fi'ri'iice'  43  are'  shown.  The'  ae  tual  isothe'rms  of  liotli  sc'a  and 
air  will  vary  throughout  I'ach  month  and  (he'ir  monthly  means  will  vary  from  year  to  yi'ar:  houe'vi'r. 
comtiaring  Figure's  23  to  28  with  Figures  29  to  .34  doe's  sugge'.st  that  the'se'  monthly  means  takc'ii  m 
combination  elo  (iri'diet  the'  ari'as  whe'ri'  icing  is  to  In'  ('X(n'cle'd  during  the'  month. 
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In  Nov(“mi)er  (Figs.  23  &  29)  no  icing  has  lu'cn  roportcci  soutli  of  4K“N,  and  all  instani-cs 
have  been  close  to  tile  Labrador  or  Newfoundland  coasts,  observations  which  accord  with  the  sea  and 
air  temperature  criteria.  That  icing  has  not  lieen  reporti-d  more  northerly  than  about  ,57  \  during 
November  is  doubtless  due  to  the  sparseness  of  Canadian  fi.shing  activity  in  these  waters;  nevertheless 
the  potential  for  serious  icing  in  the  Davis  Strait  area  exists.  Later  in  the  winter  the  southerly  advain-e 
of  sea  ice  tends  to  set  a  northern  limit  to  significant  icing  since  the  ice  cover  inhibits  wave  ai  lioii  and 
hence  spray.  Light  to  moderate  icing,  however,  may  still  be  exiierienced  wherever  sufficient  open  water 
is  encountered.  Considerable  variability  in  ice  coverage  from  year  to  year  occurs. 

In  December  (Figs.  21  is  30)  the  possibility  of  encountering  icing  moves  further  south  with 
the  advance  of  the  cold  water  and  the  lowering  air  temperatures.  Probabilities  are  still  low  in  most 
areas,  and  no  encounters  have  been  reported  in  ICN.AF  Divisions  3.M.  3N  and  4V's.  The  singli’  report  of 
icing  in  Division  3-0  was  on  the  boundary  with  3L  and  3P,  so  that  the  13C  probability  indicatc'd  may 
be  misleadingly  high. 

January  (Figs.  25  «.<:  31)  shows  considerable  incri'ase  in  icing  probabilitii's,  with  tht'  cold  water 
now  covering  all  of  the  Canailian  continental  shelf.  The  pattern  for  Hebruary  (Figs.  26  &  32  I  differs 
little  from  that  of  January  except  for  the  advance  of  sea  ic(‘  off  the  Labrador  Newfoundland  coast  and 
in  the  Ciulf  of  St.  Lawrence. 

•March  (  Figs.  27  sN  33)  sees  a  general  lessening  of  icing  probabilitii’s  as  the  air  temperature 
begins  to  moderate,  while  only  isolated  ca.ses  of  icing  have  been  reported  in  .April  (  Figs.  2iS  A;  34  ), 

7.3  .Meteorological  Conditions  Causing  Icing 

For  the  most  part  the  icing  reports  revealeil  that  icing  occurrences  followed  closely  upon  the 
passage  of  welUleveloped  low  pressure  systems  through  the  area  (Ref.  5).  These  systems  bring  in  their 
wake  cold  .Arctic  air  down  over  the  fishing  grounds  and  the  reports  indicate  that  they  induence  thi'  area 
for  only  a  matter  of  a  day  or  two  at  ;i  time.  Figures  35  and  36  illustrate  such  systems.  Occasionally, 
however,  an  extensive  outbreak  of  cold  air  will  establish  itself  over  the  eastern  part  of  the  continent  for 
a  more  extended  period,  and  then  icing  conditions  m;iy  persist  for  a  wi'ek  or  two  but  with  occasional 
lulls.  Systems  responsible  for  a  few  of  the  more  .severe  icing  conditions  in  recent  years  are  illustrated  in 
Figures  37  to  40. 

Oni'  exteniled  period  of  icing  was  from  26  January  1972  to  10  February  1972.  and  Figure  37 
shows  the  weather  pattern  iirevailing  at  1200  (LM.'l'.  on  27  January  1972.  Some  icing  was  reported  each 
day  throughout  this  periotl  with  severe  icing  being  experienced  on  26  to  2S  January.  2  Februar\’.  6  Feb¬ 
ruary  and  9  February,  so  that  many  ves.si'ls  took  precautionary  jiction  such  as  dodging  into  wind  or 
sheltering  in  the  lee  of  land.  During  these  severe  periods  air  temperatures  as  low  as  IS^C  and  winds  of 
up  to  force  11  (60  knots)  were  reported,  while  the  sea  temperature  where  many  of  the  reporting  vessels 
were  operating  was  in  the  range  0  to  3°('.  Such  i-onditions  can  result  in  very  severe  icing,  and  indeed 
this  is  home  out  by  many  of  the  reports,  since  icing  rates  in  excess  of  15  cm ''24  hr  wi're  evidenced  on 
.several  occasions.  One  stern  trawler  of  630  gro.ss  tons  reported  an  estimati'd  accumulation  of  200  tons 
of  ice  -  this  was  clearly  a  gro.ss  overi'stimate,  but  ice  thicknesses  of  u|i  to  20  inches  (51  cm  )  occurred 
despite  several  de-icings  over  a  period  of  40  hours;  the  stability  of  the  vessel  was  reported  to  have  been 
noticeably  affected.  .Another  stern  traw  ler  of  620  gross  tons  landed  at  Halifax  in  an  iced-up  condition; 
ice  from  several  locations  was  removed  and  weighed,  and  basixl  on  these  weights  and  the  thickness  of 
ice  at  various  positions  a  total  ice  weight  of  30.81  tons  was  arrived  at.  having  a  vertical  l  entre  of  gravity 
of  9.32  feet  (2.84  m)  above  deck  at  midships.  The  ice  allowance  for  this  vessel  in  its  stability  calcula¬ 
tions  is  30  tons  with  a  centroid  9.95  fei’t  (3.03  m)  above  the  level  of  the  deck  amidships.  The  greatest 
ice  thickness  on  this  vessel  was  12  inches  on  the  forward  house  front  (Fig.  41 ).  .A  government  survey 
ship  of  65  m  overall  length  sailing  from  .Newfoundlan.l  to  Halifax  in  the  same  ()eriod  picked  up  an  esti¬ 
mated  150  tons  of  ice  and  reported  an  increased  roll  pi-riod;  Figure  42  shows  some  of  this  ice  build-up 
in  the  forecastle  area. 

These  occurrences  suggest  that,  compared  with  tne  IMCO  recommendations,  the  more  severe 
ice  allowance  imposixl  by  Canada  (Ref.  44)  on  vessels  sailing  in  northern  or  eastern  ('anadian  waters  is 
not  unrealistic  for  thi'se  areas. 
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Another  extended  (jeriod  of  ieintt  was  the  12  days  from  9  February  1979  to  20  February 
1979,  with  pruticularly  severe  conditions  prevailing  on  14  and  15  February.  The'  synoptic  situation  at 
0000  G.M.T.  on  14  February  is  shown  in  Figure  39  while  that  at  1800  G.M.T.  on  15  February  is  shown 
in  Figure  40.  'I'he  more  southerly  low,  centred  at  about  40°N  55°W  at  0000  G..M.T.  on  the  14th,  deep¬ 
ened  and  moved  around  the  area  of  low  pressure  centred  just  off  the  Labrador  coast.  By  1800  G..M.'r. 
on  the  15th,  the  latter  had  dissipated  off  the  northerly  tip  of  Newfoundland  while  the  former  had 
reached  its  greatest  intensity  and  was  centred  just  south  of  Greenland.  Throughout  this  period,  cold 
air  streamed  across  Labrador  and  the  sea  areas  off  Newfoundland  and  Nova  Scotia. 

7.4  Relation  Between  Reported  Severity  and  Temperature  and  Wind  Conditions 

To  demonstrate  the  effect  of  air  temperature  and  the  relative  wind  speed  on  the  reported 
icing  severity.  Figures  43  and  44  have  been  constructed  using  two  winters’  reports  {1971-72  and  1975- 
76),  for  vessels  steaming  at  5  knots  or  greater  (Fig.  43)  and  for  vessels  making  less  than  5  knots  (Fig.  44  ). 
i.e.  w'hen  riding  out  the  storm  (dodging)  or  when  fishing.  Although  best  fit  hyperbolae  have  been  draun 
for  each  severity  level,  considerable  scatter  of  the  data  points  exists  and  no  clear-cut  transition  from  one 
severity  to  another  emerges;  nevertheless  the  trend  is  apparent.  The  scatter  may  be  attrituited  to  the 
various  factors  enumerated  in  Section  2  above  as  well  as  to  the  subjective  nature  of  the  reported  icing 
severities. 


Figures  43  and  44  may  therefore  be  regarded  as  no  more  than  rough  rules  of  thumb  in  pre¬ 
dicting  icing  seventy  if  only  the  air  temperature  and  wind  velocity  are  known.  Figure  43  .sugge.sl.s  ilia! 
when  steaming,  heavy  icing  may  be  expected  if  the  jiroduct  of  the  relative  wind  speed  (in  knots)  and 
the  air  temperature  (in  °C  below  0°C)  is  greater  than  about  400,  while  when  dodging  or  fishing  (Fig.  44  | 
the  criterion  for  heavy  icing  increases  to  a  product  in  excess  of  about  500.  Both  figure.s  suggest  that 
with  products  of  relative  wind  speed  and  air  temperature  less  than  200,  the  icing  severity  is  of  little  or 
no  significance. 

In  the  following  section,  an  attempt  is  made  to  quantitatively  predict  the  rate  of  ice  build-up 
on  a  vessel  using  an  analytical  approach. 


8.0  PREDICTION  OF  RATE  OF  ICING 

Certain  empirical  approaches  have  been  devised  to  determine  more  or  less  (juantitalively  the 
severity  of  icing  from  the  prevailing  sea  and  atmospheric  conditions;  perhaps  the  best  known  of  these 
is  due  to  Mertins  (Ref.  45).  In  his  approach.  Merlins,  in  addition  to  air  temirerature  and  winil  force, 
introduced  the  effect  of  sea-water  temperature  and  presented  his  results  in  tiie  form  of  four  charts 
(Fig.  45). 


More  recently,  Comiskey  (Ref.  46)  has  taken  Mertins’  four  charts  and  made  them  more  con¬ 
venient  to  use  by  combining  them  into  one  chart,  as  .shown  in  Figure  46.* 

.A  quantitative  rate  of  icing  may  be  approximated  from  both  the  .Mertins  and  Comiskey  charts 
by  Interpolation  using  the  range  of  icing  rates  given  for  each  icing  category. 

Applying  some  of  the  data  refiorted  by  means  of  our  icing  (pK'stionnairi's,  the  .Mertins/ 
Comiskey  charts  give  good  i-orrelation  in  those  ca.ses  when  th<‘  ve.s.sels  art"  steaming  (bt'tween  about 
5  knots  and  12  knots)  but  tend  to  underestimate  the  icing  rates  tlerived  from  tlie  re|iorts  liy  a  fadtir  of 
about  0.7.  In  cases  when  the  vessels  are  at  low  speed  (dodging  or  fishing)  correlation  is  not  so  good  and 
gt'iit'rally  the  charts  rather  seriously  overestimate  the  icing  rate.  Wind  spet'ds  relativt'  to  tin'  vt'sst'ls  were 
u-setl  rather  than  absolute  wind  speeds,  as  these  gave  somewhat  imiiroved  results. 
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Kachurin  et  al  (Ref.  47)  havo  taken  an  analytical  approach  and  calculate  the  rale  of  icing  of 
a  cylindrical  rod  under  conditions  similar  to  the  spray  icing  of  small  fishing  boats,  and  relate  this  to  th(' 
ice  build-up  rate  observed  on  actual  ships.  This  approach  apjiears  to  be  useful,  but  as  formulated  by 
Kachurin  tends  to  be  rather  complex  inasmuch  as  he  introduces  time  dependent  effects. 

Such  sophistication  seems  to  be  rather  unjustified  in  view  of  the  many  uncertainties  in  meas- 
uting  or  estimating  the  various  atmospheric  and  ocean  parameters  as  well  as  a  representative  rate  of 
icing  for  the  vessel.  Therefore  in  Reference  42  an  expression  was  developed  for  estimating  the  rale  of 
icing  on  a  ship,  wherein  only  the  primary  heat  transfer  terms  were  considered,  and  where  in  addition 
the  ice  formation  process  was  considered  to  be  a  continuous,  quasi-steady  state  process,  notwithstanding 
that  the  water  arrives  at  the  icing  surfaces  in  a  rather  discontinuous  manner. 

In  Appendix  B  of  this  report  an  amended  version  of  the  analysis  given  in  Reference  42  is 
presented  wherein  the  drop  cooling  equation  is  modifit'd  to  include  the  effect  of  evaporative  cooling, 
as  well  as  certain  other  minor  modifications. 

The  final  expression  for  the  rate  of  icing,  N,*,  derived  in  Appendix  B  is: 


N: 


*  =  0.00425  (t,  -  tj)  +  0.0369  I (t„  tj -t- 17.3  (e„  0.9e,)|mm/h 


where 


H.. 

Vu 


ta 


wave  height,  m 
relative  wind  speed,  kt 

equilibrium  temperature  of  the  icing  surface,  °C 
temperature  of  impinging  water  drops,  °C 
air  temperature,  °C 

saturation  vapour  pressure  of  moist  air  at  temperature  t^,,  t,  rt'Sjtectively.  kPa 


In  order  to  test  this  model  of  icing  under  conditions  of  sea  si)ray,  the  icing  rate,s  as  calculated 
by  this  equation  using  the  conditions  applicable  to  39  reported  icing  indicents  have  been  compared  with 
the  mean  icing  rates  derived  from  these  reports.  This  comparison  is  shown  in  Fjgiire  47;  howevf'r,  before 
discussing  this  comparison,  some  discussion  of  the  data  on  which  it  is  ba.sed  is  necessary. 

8.1  Discussion  of  Observed  Data 

The  rfclevant  observed  and  derived  data  for  39  reported  icing  incidents  in  the  two  seasons 
1977-78  and  1978-79  are  presented  as  Table  VIII.  Only  data  from  fishing  trawlers  are  used;  those  from 
freighters,  tankers  and  government  vessels  of  vairious  types  have  not  been  included  because  of  the  dif¬ 
ferent  size  and  lines  of  these  vessels,  and  until  some  means  of  incorporating  these  two  variables  into 
the  analysis  is  developed,  each  class  of  vessel  must  be  treated  separately.  Even  among  the  fishing  trawl¬ 
ers  considerable  size  and  design  variation  exists  which  undoubtedly  results  in  part  in  some  of  the  scatter 
of  the  data. 

Scatter  may  also  be  caused  by  numerous  undertainties  in  the  data,  as  the  following  discussion 

indicates. 

The  reported  wind  speed  and  direction  are  combined  with  the  vessel  speed  and  heading  to 
derive  the  wind  speed  and  direction  relative  to  the  ship.  The  relative  wind  speed  is  entered  as  column  3 
of  Table  VIII.  Data  sets  (entries)  Nos.  1  to  32  are  for  relative  wind  directions  of  0  ±  45°,  while  entries 
33  to  39  are  for  those  cases  of  beam  or  quartering  winds  (i.e.  0  >  ±45°).  Some  uncertainty  exists  iti 
interpreting  the  reported  data  since  in  many  cases  the  reported  wind  velocity  may  in  fact  l)e  the  relative 
velocity  rather  than  the  true  velocity,  and  it  is  not  always  made  clear  whether  the  ship  heading  is  mag¬ 
netic  or  true. 
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The  air  temperature  may  be  expected  to  be  fairly  reliable,  assuming  that  the  thermometer 
used  is  in  calibration  and  that  the  thermometer  bulb  is  not  exposed  to  the  water  spray.  Some  vessels 
are  evidently  not  equipped  with  thermometers,  and  so  data  from  them  could  not  be  used  except  where 
temperatures  could  be  inferred  from  other  vessels  operating  in  the  same  area  at  the  same  time. 

The  sea  water  temperature  entered  in  column  5  of  'Fable  Vlll  is  derived,  knowing  the  vessel 
location,  from  weekly  sea-surface  temperature  charts  compiled  by  the  Canadian  Forces  Meteorology/ 
Oceanography  Centre,  Halifax. 

The  wave  height  is  not  always  quantitatively  reported  and  in  most  instances  has  been  inferred 
from  a  qualitative  description  of  the  sea  state  using  the  following  scale  derived  from  the  State  of  Sea 
table  in  the  Trawlermen’s  Handbook  (Ref.  48). 


Description 

Inferred  Wave 
Height  —  m 

slight 

0.6  -  1.2 

moderate 

1.2  -  2.5 

rough 

2.5  -  4 

very  rough 

4  -  6 

high 

6  -  9 

very  high 

9  -  14 

If,  as  stated  in  Appendix  B,  the  liquid  water  content  of  the  spray  is  primarily  a  function  of 
wave  height,  then  this  is  an  important  parameter  in  the  analysi.s,  but  one  which  is  sultject  to  consideralile 
uncertainty  because  of  the  necessity  of  interpreting  the  qualitative  descriptions  given  in  the  reports. 
Relationships  relating  wave  height  to  wind  speed  exist,  but  these  usually  apply  to  winds  of  unlimited 
duration  and  fetch,  conditions  which  do  not  apply  to  the  majority  of  the  fishing  trawler  icing  events 
liecause  of  their  proximity  to  land  or  to  an  ice  cover. 

'Fhe  salinity  of  the  sea  water  in  the  location  of  an  icing  report  has  been  derived  from  the 
.Mean  Sea  Surface  Salinity  charts  given  in  the  Oceanographic  Atlas  (Ref.  37).  Bc'cause  changes  in  salinity 
within  the  range  to  he  ('xpected  in  this  area  (i.e.  30%..  <  s  <  35/fo)  do  not  greatly  affect  the  icing  rati', 
secular  variations  from  the  salinity  values  given  in  the  charts  are  of  little  significance. 

Perhaps  the  greatest  uncertainty  arises  in  the  derivation  of  icing  rate  (column  8).  A  mean  or 
ri'iiresentative  ice  thickness  was  assessed  from  the  thicknes,ses  reported  for  the  various  on-board  locations 
indicated  in  the  Icing  Report  (Fig.  19).  In  a.ssessing  the  mean  thickness,  somewhat  less  weight  was  ac¬ 
corded  the  ice  thickness  on  the  forward  rails  and  stay,  since  these  thicknesses  were  not  nece.ssarily  rep¬ 
resentative  of  the  vessel  as  a  whole.  The  reporting  of  duration  of  icing  occasionally  left  .some  uncertainty, 
particularly  when  icing  was  intermittent,  in  which  case  what  should  have  lieen  reported  as  seiiarate  (‘vents 
were  .sometimes  reported  as  one  extended  event.  Sometimes,  too,  it  was  not  made  clear  whether  some 
de-icing  action  that  could  have  affected  the  reported  ice  thicknesses,  had  l)een  taken  during  the  icing 
encounter. 

Other  undertainties  arise  because  of  changing  conditions  during  the  encounter.  Often  such 
changes  (in  wind,  temperature,  etc.)  are  not  indicated.  Olearly,  the  longer  the  reported  encounter,  the 
more  uncertainties  there  are,  so  that  no  encounter  of  duration  greater  than  24  hours  has  lieen  used  in 
Talile  VIII. 

8.2  Discussion  of  Icing  Rate  Results 

Kmploying  Kquation  (B.31c)  of  Afipcndix  Band  the*  iterative  techniques  described  (Section 
B.1.7),  the  data  of  columns  2  to  7  of  'Fable  VIII  were  u.sed  to  obtain  calculated  rates  of  icing  (column  9 
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of  Table  VIII).  These  are  compared  in  Figure  47  with  the  rates  of  icing  as  derived  from  the  icing 
reports  (column  8  of  Table  VIII).  Also  shown  in  the  figure  is  the  line  passing  through  the  origin  that 
gives  the  best  least  squares  fit  between  the  two  sets  of  values.  The  correlation  coefficient  is  r  =  0.55. 

That  there  is  so  much  scatter  in  Figure  47  is  perhaps  not  unexi)ected  in  view  of  the  numerous 
uncertainties  in  the  raw  data  discussed  in  sub-section  8.1.  It  is  encouraging  in  fact  that  the  calculated 
rate  is  generally  within  an  order  of  magnitude  of  the  assessed  rat»‘.  which  gives  a  certain  degree  of  con¬ 
fidence  in  the  calculations;  however,  a  more  complete  test  of  the  theory  cannot  be  made  until  more 
consistent  and  reliable  data  is  available.  That  some  of  the  present  tlata  is  “bad”  may  be  seen  by  com¬ 
paring  some  of  the  entries  in  Table  VIll,  c.f.  entries  1  and  15,  25  and  33,  7  and  15,  10  and  21,  etc. 
where  lower  reported  icing  rates  correspond  to  clearly  more*  severe  conditii)ns. 

Figure  47  does  not  suggest  that  a  beam  wind  results  in  a  significantly  different  rate  of  icing 
than  a  bow  wind,  nor  can  any  consistent  effect  of  ve'ssel  length  he  observc'd. 

The  rate  of  icing  Equation  (B.31c)  may  be  used  to  show  the  relative  effects  of  the  various 
parameters  contributing  to  the  heat  balance.  This  is  illustrated  in  Figuri*  48  where  one  parameter  of 
a  set  of  standard  conditions  has  been  varied  in  turn  throughout  its  likely  range.  The  standard  conditions 
chosen  were: 


Air  temperature 

T_,  = 

lO'^C 

Sea-surface  temperature 

o°c 

Relative  wind  s|)eed 

V  = 

40  kt 

Wave  height 

= 

3  m 

Salinity 

s 

32'ro 

Relative  humidity 

Rn  = 

90^-; 

Figure  48  indicates  that  the  most  significant  variables  in  determining  the  icing  .severity  are 
the  air  temperature  and  the  relative  wind  speed.  The  sea-surface  temperature,  the  wave  height  and 
relative  humidity  are  seen  to  be  of  more  moderate  significance,  while  the  salinity  of  the  sea  water  is 
shown  to  have  the  least  effect  on  rate  of  icing  of  these  six  variables. 

The  results  demonstrated  in  Figure  48  therefore  suggest  that,  in  view  of  the  quality  of  the 
available  data,  a  simplification  may  be  made  to  the  Rate  of  Icing  equation  by  assuming  a  constant 
■salinity  and  a  constant  wave  height  (i.e.  a  constant  water  content,  say  0.0005  kg/m’ ).  Further  simplifica¬ 
tion  may  lie  made  by  assuming  that  the  icing  surface  temperature  is  identical  to  the  freezing  tempera¬ 
ture  of  the  sea  water.  If  then  the  icing  surface  temperature  is  taken  as  1.8°(',  there  results  the  sim¬ 
plified  equation: 


Ni*  =  0.013  (  t, I  -  1.8)V  -t-  0.0369V‘"*  |(-t.,  1.8)  17.3  (0..5354  0.9  e , ) 


in  which  only  T^,  has  to  be  determined  from  Equation  (B.14)  and  e^,  by  Equation  (B.28). 

The  application  of  this  equation  to  the  data  of  Table  VIII  results  in  a  plot  not  significantly 
differi'nt  from  Figure  47  (i.e.  a  correlation  coefficient  of  0.57  and  the  best  fit  lim*  through  the  origin 
having  a  slope  of  0.74). 


9.0  CLOSURE 

This  report  has  collated  and  reviewed  the  work  in  the  programme  on  the  icing  of  fishing 
trawlers  that  the  Low  'I’emperature  Laboratory  has  done  over  the  eleven-year  pi*riod  1969  to  1979. 
Because  of  budgi'tary  restraints  this  work  has  been  very  low  key  in  nature*. 


Ml 
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Initially,  some  simple  icing  tc'sts  (Ref.  17)  on  cylinders  of  various  diameters  were  made  in  the 
icing  wind  tunnel  to  demonstrate  the  effects  of  air  temperature  and  cylinder  size  on  the  resulting  ice 
formations.  The  main  conclusion  from  these  tests  was  that  the  icing  efficiency  of  a  cylinder  increases 
both  with  decreasing  diameter  and  with  decreasing  tem{)erature,  a  result  that  demonstrates  the  need  to 
reduce  small  scale  components  such  as  spars,  rigging,  etc.  to  a  minimum,  and  where  they  cannot  l)e 
eliminated  cdtogether,  wherever  possible  to  let  one  larger  diameter  component  serve  the  purpose  of 
several  smaller  ones. 

Among  a  number  of  simple  methods  tested  as  possible  means  of  reducing  the  icing  hazard 
(Ref.  27),  no  passive  method  was  found  to  significantly  ea.se  the  removal  of  ice  from  a  surface.  However, 
an  inflatable  rubber  de-icing  lilanket  was  found  to  be  remarkably  effective  in  removing  ice  from  both 
cylindrical  and  flat  surfaces. 

Icing  report  forms  have  been  distributed  since  1969  to  vessels  operating  out  of  eastern  Cana¬ 
dian  ports.  The  data  olitained  by  this  means  has  helped  to  establish  the  type  of  weather  systems 
responsible  for  icing  conditions,  the  geographical  extent  of  the*  occurrence  of  icing,  and  statistics  on 
the  severity  of  icing  encounters  in  the  study  area. 

The  geographical  extent  varies  from  month  to  month,  l)ut  at  its  most  e.xtensive  covers  the 
whole  of  the  continental  shelf  off  Eastern  Canada  as  well  as  the  Labrador  Sea.  Tlie  type  of  weather 
system  mostly  responsilde  for  severe  icing  in  these  areas  is  a  well-develop(>ii  low  pressure  system,  in  the 
wake  of  which  cold  Arctic  air  is  brought  down  over  the  fishing  grounds  by  vigorous  northerly  or  nortli- 
westerly  winds. 

The  reports  show  that  during  January  and  February  icing  conditions  of  light  or  greater 
intensity  occur  on  the  average  one  day  out  of  two  somewhere  in  the  study  area,  while  in  bad  years  this 
may  increase  to  nearly  three  days  in  four.  In  terms  of  the  amount  of  ice  accumulated,  it  is  shown  that 
on  approximately  8  percent  of  voyages  during  the  winter  season,  ice  accumulations  in  excess  of  the 
allowance  recommended  by  IMCO  in  calculating  the  vessel’s  staljility  were  experienced. 

Finally,  an  analytical  approach  lo  the  prediction  of  the  rate  of  ic('  l)uild-ui)  on  a  fishing 
trawler  is  presented.  This  analysis  (with  possil)le  minor  modifications)  has  application  also  to  the  icing 
of  marine  structures. 
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TABLE  I 


(Reproduced  from  Ref.  38) 


Beaufort 
Force  of 
Wind 


Theoretical 
Max.  Wave 
Height  (ft), 
Unlimited 
Duration 
and  Fetch 


Duration  of  Winds  ( hours) 
with  Unlimited  Fetch,  to 
Produce  Percent  of  Max. 
Wave  Height  Indicated 


Fetch  (n.  miles),  with 
Unlimited  Duration  of 
Blow,  to  Produce  Percent 
of  Max.  Wave  Height 
Indicated 


50% 

90^^ 

50% 

7^ 

90% 

3 

2 

1.5 

5 

8 

3 

13 

25 

5 

8 

3.5 

8 

12 

10 

30 

60 

7 

20 

5.5 

12 

21 

22 

75 

150 

9 

40 

7 

16 

25 

55 

150 

280 

11 

70 

9 

19 

32 

85 

200 

150 

PAufc  bamk-mot  rij-ett 
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TABLE  II 

ICING  QUESTIONNAIRE  RESPONSE  BY  YEAR 


Number  of  Vessels  Responding  and  Number  of  Reports 
(including  Null  Reports)  Received^®^ 


Year 

Fishing 

Vessels 

Government 

Vessels’ 

Freighters, 
Tankers,  etc. 

Ferry 

Services 

A 

B 

A 

B 

A 

B 

A  B 

1968-69 

2 

2 

3 

3 

1 

1 

- 

1969-70 

13 

28 

— 

- 

1 

1 

- 

1970-71 

29 

89 

- 

-■ 

- 

- 

1971-71 

24 

89 

1 

1 

3 

13 

1  3 

1972-73^ 

■- 

- 

- 

-■ 

1  5 

1973-74 

29 

119 

-- 

- 

1 

2 

1  1 

1974-75’ 

5 

8 

15 

43 

8 

8 

1  40 

1975-76 

38 

172 

9 

29 

1 

2 

2  52 

1976-77 

35 

124 

10 

28 

3 

5 

4  70 

1977-78 

48 

182 

17 

76 

1 

2 

1  4 

1978-79 

22 

73 

13 

20 

'I’otals 

886 

200 

34 

175 

Notes: 

1  ( lovernmetit  vessi-l.s  m(  lii<le  nc  l)reakers.  survey  vessels,  supply  vessels,  fisherv’  patrol  vessels,  eti 

2  Itt  year  l!>72  7.'l  reports  from  N'ewfountlland  were  lost. 

.1  111  year  1  !t7  I  of)  Newfoimdland  fishing  flei-ts  wen-  on  strike 


TABLE  III 


ICINC;  SEVERITY  STATISTICS  FOR  THE  WHOLE  STUOY  AREA 


No.  of  Days  Icing  Severity  Reported  as: 


Light 
or  Oreater 

Moderate 
or  (Ireater 

Heavy 

December 

7-year  total 

32 

26 

0 

worst  year 

10 

10 

3 

yearly  avt-ragt' 

l.G 

3.7 

1,3 

standard  deviation 

3.5 

3.2 

1.0 

daily  ex()eetation 

M.7'; 

12.0'; 

1,1'- 

January 

7-year  tcital 

OH 

70 

21 

worst  year 

23 

15 

6 

yearly  average 

1  1.0 

10.0 

3.0 

standard  tlc'viation 

1.5 

3.0 

1.6 

daily  expe(.'tation 

15.2'; 

32.2"' 

0.7'; 

February 

7-year  total 

101 

70 

32 

worst  year 

20 

1(1 

I 

yearly  averugt- 

1  1.1 

10,0 

1.6 

standard  deviation 

3.1 

3.7 

2.1 

daily  exfiectation 

51.0'; 

35. }'; 

16.2'- 

March 

7 -year  total 

71 

17 

11 

worst  year 

27 

17 

; ) 

yearly  average 

10,1 

<>.7 

2.0 

standard  devitition 

H  !> 

5.S 

l.!( 

daily  expectation 

32.7'; 

21 .7'- 

6.5'- 

*  *  * 

*  ♦  * 

1  l)e<'eml>er  -  .11  March 

7  year  total 

302 

2 1 .3 

76 

worst  year 

72 

50 

15 

yearly  average 

13  1 

30  1 

10.0 

standard  deviation 

13  0 

10  5 

3.1 

dally  expectation 

35  tv- 

25.1', 

0.0'- 
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TABLE  IV 

ICING  SEVERITY  STATISTICS  FOR  ICNAF  DIVISION  3L 


No.  of  Days  Icing  Severity  Reported  a.s: 


Light 
or  Greater 

Moderate 
or  Greater 

Heavy 

December 

7 -year  total 

7 

5 

5 

worst  year 

3 

3 

3 

yearly  average 

1.0 

0.7 

0.7 

standard  deviation 

1.1 

1.2 

1.2 

daily  expectation 

S.2'i 

2.3'; 

2.3'; 

January 

7 -year  total 

33 

24 

9 

worst  year 

11 

8 

3 

yearly  average 

1.7 

3.4 

1.3 

standard  deviation 

2.0 

2.1 

1.3 

daily  ('xpt'ctation 

15.2'; 

11. r; 

i.r; 

February 

7 -year  total 

30 

11 

2 

worst  year 

8 

4 

1 

yearly  average 

1.3 

1.(1 

0.3 

standard  deviation 

1.9 

1.4 

0.5 

daily  expectation 

15.2'; 

5.6'; 

1 .0'; 

March 

7 -year  total 

27 

21 

4 

worst  year 

12 

8 

2 

yearly  average 

3.9 

3.0 

0.6 

standard  deviation 

3.9 

2.9 

0.7 

daily  expectation 

12.4'r 

9.7''; 

1.8'; 

*  *  ♦ 

♦  * 

*  ♦  * 

1  December  —  31  March 

7 -year  total 

97 

61 

20 

worst  year 

34 

20 

5 

yearly  average 

13.9 

8.7 

2.9 

standard  deviation 

8.8 

5.3 

1.8 

daily  expectation 

11.4'; 

7.2'. 

2.4'; 
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TABLK  V 

ICING  SKVKRITY  STATISTICS  FOR  ICNAF  DIVISION  3P 


No.  of  Days  Icing  Severity  Reported 

as; 

Light 
or  Greater 

.Moderate 
or  Greater 

Heavy 

Dofenibcr 

7-year  total 

4 

3 

0 

worst  year 

2 

2 

0 

yearly  averaj^e 

0.6 

0.4 

0 

stanciard  deviation 

0.7 

0.7 

0 

daily  expectation 

1.8'; 

1  A'< 

O'; 

January 

7-year  total 

53 

38 

10 

worst  year 

11 

8 

3 

yearly  average 

7.6 

5.4 

1.4 

standard  deviation 

3.0 

2.7 

1.2 

dialy  expectation 

24.4'; 

17,5'; 

4.6'; 

February 

7-year  total 

51) 

41 

18 

worst  year 

16 

12 

4 

yearly  average 

8.4 

5.9 

2.6 

standard  deviation 

3.6 

3,1 

1.3 

daily  expei-tation 

20.8'; 

20.7'; 

9.1'; 

March 

7-year  total 

38 

23 

4 

worst  year 

14 

7 

3 

yearly  average 

5.4 

3.3 

0,6 

standard  deviation 

5.2 

2.8 

1.0 

daily  expectation 

1  7.5''r 

10.6'; 

1 .8'V 

*  ♦ 

*  *  * 

♦  *  * 

1  December  —  31  March 

7-year  total 

1 54 

105 

32 

worst  year 

37 

27 

7 

yearly  average 

22.0 

15.0 

4.5 

standard  deviation 

8,3 

6.3 

1.7 

daily  ('xpectation 

I8,r; 

12.4'. 

3.8'; 

SUMMARY  OF  ICING  ENCOUM  KRS  IN  I'ERMS  OF  IMCO  STANDARD  ICINX;  ALLOWANCE 
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TABLE  VII 

NUMBER  (AND  PERCENTAGE)  OF  YEARS  FOR  WHICH  AT  LEAST  ONE  OCCURRENCE  OF 
ICING  OF  SEVERITY  GREATER  THAN  SLIGHT  HAS  BEEN  REPORTED  FOR  THE 
DESIGNATED  ICNAF  DIVISION  IN  EACH  OF  THE  WINTER  MONTHS 


ICNAF 

Division 

Nov. 

Dec. 

Jan. 

Feb. 

Mar. 

Apr. 

IF 

1  (13%) 

2G 

1  (13%) 

2H 

1  (13%) 

2  (25%) 

1  (13%) 

2J 

3  (38%) 

3  (38%) 

3  (38%) 

1  (13%) 

3K 

2  (25%) 

1  (13%) 

4  (50%) 

4  (50%) 

2  (25%) 

3L 

1  (13%) 

2  (25%) 

7  (88%) 

7  (88%) 

4  (50%) 

3M 

3N 

1  (13%) 

2  (25%) 

30 

1  (13%) 

2  (25%) 

5  (63%) 

2  (25%) 

3P 

3  (38%) 

8  (100%) 

8 (100%) 

6  (75‘/ ) 

2  (25%) 

4R 

1  (13%) 

2  (25%) 

5  (63%) 

4  (75%) 

1  (13%) 

2  (25%) 

4S 

3  (38%) 

4  (50%) 

1  (13%) 

1  (13%) 

4T 

2  (25%) 

2  (25%) 

3  (38%) 

1  (13%) 

4Vn* 

4  (50%) 

8  (100%) 

6  (75%) 

3  (38%) 

2  (25%) 

4Vs* 

4W 

4  (50%) 

4  (50%) 

7  (88%) 

2  (25%) 

4X 

3  (38%) 

5  (63%) 

4  (50%) 

3  (38%) 

5Y 

1  (13%) 

1  (13%) 

2  (25%) 

1  (13%) 

Bfised  on  8  years  of  data 


* 


The  boundary  between  4Vn  and  4Vs  is  here  taken  to  be  the  parallel  44°  10'  N  instead  of  45°  40’  N. 
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TABLE  VIII 

DATA  USABLE  FOR  ICING  RATE  STUDIES 
FROM  SEASONS  1977-78  AND  1978-79 


(1) 

Data 

Set 

No. 

(2) 

Vessel 

Length 

m 

(3) 

Rel.  Wind 
Velocity 
kt 

(4) 

Air 

Temp. 

°C 

(5) 

Sea-Surf. 

Temp. 

°C 

(6) 

Wave 

Height 

m 

(7) 

Sea 

Salinity 

%o 

(8) 

Icing 

Rate 

mm/h 

(9) 

Calculated 
Icing  Rate 
(Eq.  B.3/C) 
mm/hr 

1 

51.5 

55 

-6 

+1 

5* 

32.25 

8.5 

5.1 

2 

51 

35 

-3.5+ 

-0.5 

1.2 

31 

1.6 

2.1 

3 

51 

35 

-6  t 

+4 

3* 

32.35 

0.79 

3.3 

4 

51 

48 

-3 

-0.5 

4.5* 

34 

4.2 

1.6 

5 

51 

32 

-5 

-1 

1 

33.75 

1.3 

3.0 

6 

51 

56 

-4 

-2 

3* 

33.25 

0.75 

4.2 

7 

51 

36 

-7  t 

+1 

3* 

32.25 

10.9 

5.3 

8 

51 

39 

-4 

0 

3* 

32.7 

3.2 

2.6 

9 

51 

32 

-10  t 

+2 

9 

33.75 

4.8 

7.3 

10 

47 

36 

-2 

0 

0.6* 

31.5 

4.8 

0.65 

11 

47 

25 

-5 

0 

2.4* 

34 

0.85 

2.7 

12 

47 

40 

-9  + 

0 

2.4 

32.5 

0.71 

7.9 

13 

47 

40 

-6.5 

+0.5 

3* 

32.75 

2.8 

5.2 

14 

47 

50 

-4 

-1 

3* 

32.3 

5.4 

3.5 

15 

47 

55 

-8 

+1 

3* 

32.3 

1.6 

8.2 

16 

46 

38 

-8  t 

+1.5 

2 

31 

3.8 

6.3 

17 

44 

45 

-14 

-0.5 

3* 

32.75 

12.7 

14.1 

18 

44 

35 

-5.5 

+4 

1 

32.1 

3.0 

3.3 

19 

44 

28 

-3 

+1 

2* 

32 

1.8 

1.2 

20 

44 

28 

-10 

0 

2* 

31 

7.3 

6.8 

21 

42 

37 

-10 

-0.5 

2* 

33.5 

3.6 

8.2 

22 

42 

26 

-5 

0 

1* 

31 

2.6 

2.6 

23 

42 

25 

-5 

0 

1* 

33.5 

2.1 

2.4 

24 

42 

25 

-5 

+3 

3* 

34 

1.9 

2.2 

25 

42 

25 

-6 

0 

1.5 

31 

8.0 

3.4 

26 

42 

33 

-10 

+1.5 

2 

32.5 

3.5 

7.3 

27 

42 

38 

-12 

+2 

1.5 

32.25 

6.5 

9.5 

28 

42 

30 

-10 

0 

1.2 

32 

6.9 

6.6 

29 

42 

39 

-16 

0 

1.2 

33.75 

10.6 

12.8 

30 

38 

41 

-7 

-1 

3* 

32.25 

1.6 

6.4 

31 

34 

26 

-10  + 

0 

3* 

31 

2.4 

6.8 

32 

34 

32 

-7  + 

-1 

3* 

32 

2.2 

5.4 
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TABLE  VIII  (Cont’d) 

DATA  USABLE  FOR  ICING  RATE  STUDIES 
FROM  SEASONS  1977-78  AND  1978-79 


(1) 

Data 

Set 

No. 

(2) 

Vessel 

Length 

m 

(3) 

Rel.  Wind 
Velocity 
kt 

(4) 

Air 

Temp. 

°C 

(5) 

Sea-Surf. 

Temp. 

°C 

(G) 

Wave 

Height 

m 

(7) 

Sea 

Salinity 

%o 

(8) 

Icing 

Rate 

mm/h 

(9) 

Calculated 
Icing  Rate 
(Eq.B.3/c) 
mm/hr 

33 

51 

46 

-6  t 

0 

2 

32.25 

4.0 

5.3 

34 

51 

21 

-10 

-0.5 

3* 

33.25 

3.9 

5.8 

35 

47 

29 

-5 

+0.5 

2 

33.5 

2.8 

2.9 

36 

44 

31 

-6 

0 

2* 

32.3 

1.9 

4.0 

37 

42 

35 

-6 

-0.5 

3* 

32.5 

2.1 

4.6 

38 

42 

19 

-4 

0 

2* 

32.2 

1.6 

1.7 

39 

38 

31 

-5 

+2 

3 

31 

2.5 

2.9 

Data  in  Columns  (3),  (4),(6)  and  (8)  derived  from  Ship  Icing  Reports 

t  These  air  temperatures  inferred  from  other  vessels  in  same  general  vicinity. 

*  These  wave  heights  inferred  from  sea  state  description. 


IN  DATA  SETS  1  TO  32  THE  RELATIVE  WIND  DIRECTION  IS  0°  ±45°. 
DATA  SETS  33  to  39  HAVE  RELATIVE  WIND  DIRECTION  >  45°. 


DISTANCE  BETWEEN  UPPER  AND  LOWER  TANGENT  TRAJECTORIES  IN 


JitECEiJiiJU  PAu*  BUNK -WOT  TlufELi 


FIG.  1 :  WATER  DROP  TRAJECTORIES  IN  THE  VICINITY  OF  A  CIRCULAR  CYLINDER 


FIG.  2;  ICE  ACCRETIONS  ON  15.2  CM  HORIZONTAL  CYLINDER  AFTER  1  HOUR  EXPOSURI 


CCRETIONS  ON  15  2  CM  VERTICAL  CYLINDER  AFTER  1  HOUR  EXPOSURE 

RUN  N(1  16 
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T 


6 

<j 

O' 


10-20%  LESS 
THAN  MIDSECTION 


MIOSECTION 


10-20%  GREATER 
THAN  MIDSECTION 


FIG.  4:  DIMENSIONS  OF  ICE  ON  VERTICAL  AND  HORIZONTAL  CYLINDERS 


ICING  EFFICIENCY  % 
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AIR  TEMPERATURE  ®C 


FIG.  5:  EFFECT  OF  TEMPERATURE  ON  ICING  EFFICIENCY 


FIG.  8:  TWO  PHASE  THERMOSYPHON 


G  A  ij  vj  E 

I —  OPff’ 

^  FF  vAivt 

PRESSURE  REf^uCiNG  vAly/F 

U>1  CONTHOLCiEAR 


FIG.  10:  DETAILS  OF  PNEUMATIC  DE  ICER 


TEMPERATURE:  lO  'C  ICE  THICKNESS:  UP  TO  9  cm  DE  ICER  AIR  PRESSURE.  70  kPa 


r 


(a)  AT  CONCLUSION  OF  ICING  (b)  AFTE R  TWO  I NF LATIONS  OF  DE  ICE R 


(c)  AFTER  TWO  inflations  OF  DE  ICER 


(d)  AFTER  THIRD  INFLATION 


FIG.  14:  OUTDOOR  TEST  NO.  2 


ICING  TEMPERATURE:  15°C 
DE  ICED  AT  -ig^C  AFTER  17  HOURS  COLD  SOAK 


U)  AT  CONCLUSION  OF  ICING 


(bl  AFTER  OPERATION  OF  DE  ICERS 


FIG.  15;  OUTDOOR  TEST  NO.  6 


ICING  TEMPERATURE  15"C 
DEICED  AT  12"’C 


FIG.  16:  PRINCIPLE  TRACKS  OF  INTENSE  STORMS  FOR  JANUARY, 
FEBRUARY  AND  MARCH  1963-67 


(FROM  ARCHIBALD,  Ref.  36) 


FIG.  17;  MEAN  SEA  SURFACE  TEMPERATURE  (  F)  FOR  FEBRUARY 

(FROM  OCEANOGRAPHIC  ATLAS  REF.  371 
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FIG.  18;  FREEZING  TEMPERATURE  OF  SEA  WATER  AS  A  FUNCTION  OF  SALINITY 
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SHIP  ICING  REPORT 


Pi  I  AM  MIPOHI  l)N  IMIS  K)IW  ANY  iNSIANl.fc  0^  K  t  ♦UHMAilON  i)N  VtSbtl  DUHIN'i  VOYA(A.  .  il 
Mom  IHAN  UNI  0<  tUHl’INl.  I  0»  iL-INf,  ON  VOYA(>l.  iJbl  SI  PAHA  II  I  OIW  I  QH  I  Ai  M  INI  UUNTtH 
n  »■  IMPOHIAM  lO  I  lA'H’llU  AIOHMIOHIAIH  VOYAC.I  IVIN  if  NO  I  CM .  O'OJUHb. 


'  OMPl  I  Tt  f)  f  OHM'  )l  [)  HI  MAH  I  O  TO  .  'A  hai.  mi  .  ahc.ma  'i  .m>  i  .i  . ,n  -  ,i  vn  ..  a  m,  a  t  n-  .  «.%■ 

A  ’  !•  >Tj  A  ,  HI  M  AM.  .  I  I  •  .  •  t  A  AA  .  i*j  I  AM  .  !  •  >  A  .H« 


FIG.  19:  SHIP  ICIN'  REPORT  FORM 


FIG.  22:  DESIGNATED  SEA  AREAS  OF  THE  INTERNATIONAL  COMMISSION  FOR 
THE  NORTHWEST  ATLANTIC  FISHERIES  (ICNAF) 


FIG.  23:  PERCENTAGE  OF  YEARS  FOR  WHICH  ICING  OF  SEVERITY  GREATER  THAN 
SLIGHT  HAS  BEEN  REPORTED  IN  THE  VARIOUS  ICNAF  DIVISIONS  FOR  THE 

MONTH  OF  NOVEMBER 


NUMBFR  OF  YEARS  OF  DATA:  8 


FIG.  24;  PERCENTAGE  OF  YEARS  FOR  WHICH  ICING  OF  SEVERITY  GREATER  THAN 
SLIGHT  HAS  BEEN  REPORTED  IN  THE  VARIOUS  ICNAF  DIVISIONS  FOR  THE 

MONTH  OF  DECEMBER 

Yt  ARSOF  DMA  8 


■  Gl  • 


FIG.  25:  PERCENTAGE  OF  YEARS  FOR  WHICH  ICING  OF  SEVERITY  GREATER  THAN 
SLIGHT  HAS  BEEN  REPORTED  IN  THE  VARIOUS  ICNAF  DIVISIONS  FOR  THE 

MONTH  OF  JANUARY 


NUIVlHf  H  OF  Yt  Ah’S  Of  DAI  A  8 


HIG.  26:  PERCENTAGE  OF  YEARS  FOR  WHICH  ICING  OF  SEVERITY  GREATER  THAN 
SLIGHT  HAS  BEEN  REPORTED  IN  THE  VARIOUS  ICNAF  DIVISIONS  FOR  THE 

MONTH  OF  FEBRUARY 


NUMBER  OF  YEARS  OF  DATA  8 


FIG.  27:  PERCENTAGE  OF  YEARS  FOR  WHICH  ICING  OF  SEVERITY  GREATER  THAN 
SLIGHT  HAS  BEEN  REPORTED  IN  THE  VARIOUS  ICNAF  DIVISIONS  FOR  THE 

MONTH  OF  MARCH 

NUMBER  OF  YEARS  OF  DATA:  8 


FIG.  28.  PERCEIMTAGE  OF  YEARS  FOR  WHICH  ICING  OF  SEVERITY  GREATER  THAN 
SLIGHT  HAS  BEEN  REPORTED  IN  THE  VARIOUS  ICNAF  DIVISIONS  FOR  THE 

MONTH  OF  APRIL 

NUMBER  OF  YEARS  OF  DATA:  8 


-  +4,4°C  MEAN  SEA  SURFACE  ISOTHERM 

-  10%  PROBABILITY  OF  AIR  TEMPERATURE  0°C  OR  LESS 

- APPROXIMATE  MEAN  ICE  LIMIT 

FIG.  29:  MEAN  SEA  SURFACE  TEMPERATURE  OF  40°  F  {4.4°C)  and  10%  FREQUENCY 
OF  OCCURRENCE  OF  SURFACE  AIR  TEMPERATURE  OF  O'  C  OR  LESS 
FOR  THE  MONTH  OF  NOVEMBER 


A 
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-  +4.4°C  MEAN  SEA  SURFACE  ISOTHERM 

- 10%  PROBABILITY  OF  AIR  TEMPERATURE  0°COR  LESS 

- APPROXIMATE  MEAN  ICE  LIMIT 

FIG.  30:  MEAN  SEA  SURFACE  TEMPERATURE  OF  40°F  (4.4°C)  and  10%  FREQUENCY 
OF  OCCURRENCE  OF  SURFACE  AIR  TEMPERATURE  OF  0°C  OR  LESS 
FOR  THE  MONTH  OF  DECEMBER 


+4,4  C  MEAN  SEA  SURFACE  ISOTHERM 


FIG.  31 


10%  PROBABILITY  OF  AIR  TEMPERATURE  0°C  OR  LESS 
APPROXIMATE  MEAN  ICE  LIMIT 


:  MEAN  SEA  SURFACE  TEMPERATURE  OF  40°  F  (4.4°C)  and  10%  FREQUENCY 
OF  OCCURRENCE  OF  SURFACE  AIR  TEMPERATURE  OF  0"C  OR  LESS 
FOR  THE  MONTH  OF  JANUARY 


-  +4.4°C  MEAN  SEA  SURFACE  ISOTHERM 

-  10%  PROBABILITY  OF  AIR  TEMPERATURE  0°C  OR  LESS 

- APPROXIMATE  MEAN  ICE  LIMIT 


FIG.  32;  MEAN  SEA  SURFACE  TEMPERATURE  OF  40°F  {4.4°C)  and  10%  FREQUENCY 
OF  OCCURRENCE  OF  SURFACE  AIR  TEMPERATURE  OF  0°C  OR  LESS 
FOR  THE  MONTH  OF  FEBRUARY 


W 


-  M.4°C  MEAN  SEA  SURFACE  ISOTHERM 

-  10%  PROBABILITY  OF  AIR  TEMPERATURE  0°C  OR  LESS 

- - APPROXIMATE  MEAN  ICE  LIMIT 

FIG.  33:  MEAN  SEA  SURFACE  TEMPERATURE  OF  40°F  (4.4°C)  and  10%  FREQUENCY 
OF  OCCURRENCE  OF  SURFACE  AIR  TEMPERATURE  OF  0°C  OR  LESS 
FOR  THE  MONTH  OF  MARCH 


r 


+4.4  C  MEAN  SEA  SURFACE  ISOTHERM 


10%  PROBABILITY  OF  AIR  TEMPERATURE  0°C  OR  LESS 
APPROXIMATE  MEAN  ICE  LIMIT 


FIG.  34:  MEAN  SEA  SURFACE  TEMPERATURE  OF  40°F  (4.4°C)  and  10%  FREQUENCY 
OF  OCCURRENCE  OF  SURFACE  AIR  TEMPERATURE  OF  0°C  OR  LESS 
FOR  THE  MONTH  OF  APRIL 


X  -  VESSELS  EXPERIENCING  ICING 
•  VESSEL  NOT  EXPERIENCING  ICING 

FIG.  35:  SURFACE  CHART  SHOWING  METEOROLOGICAL  SYSTEM  RESPONSIBLE  FOR 
ICING  OCCURRENCES  ON  14  FEBRUARY.  1970 


X  -  VESSELS  EKPERIENCINC  ICING 
■  -  VESSEL  NOT  EXPERIENCING  ICING 

FIG.  36:  SURFACE  CHART  SHOWING  METEOROLOGICAL  SYSTEM  RESPONSIBLE  FOR 
ICING  OCCURRENCES  ON  17  FEBRUARY,  1970 


X  VESSELS  EXPfRIFNCING  ICING 

■  VESSELS  NCn  EXPERIENCING  ICING 


SHELTERING  IN  LEE  OE  LANO 

FIG.  37  surface  WEATHER  CHART  FOR  1200  GMT,  27  JANUARY,  1972 


k 


X  VESSELS  EXPERIENCING  ICING 


■  VESSEL  IN  PACK  ICE  -  NOT  EXPERIENCING  ICING 

FIG.  38:  SURFACE  WEATHER  CHART  FOR  0600  GMT,  12  MARCH,  1974 


X  -  VESSELS  EXPERIENCING  ICING 


FIG.  39.  SURFACE  WEATHER  CHART  FOR  0000  GMT,  14  FEBRUARY,  1979 


X  VESSELS  EXPERIENCING  ICING 


FIG.  40;  SURFACE  WEATHER  CHART  FOR  1800  GMT,  15  FEBRUARY,  1979 


(a)  ICE  BUILD  UP  ON  RAILS,  ANCHOR  AND  FORESTAY 


(h)  ICE  BUILD  UP  ON  WHEELHOUSE  FRONT 


FIG,  41:  ICING  OF  STERN  TRAWLER  26  27  JANUARY  1972 


RELATIVE  WIND  SPEED  -  KNOTS 


1 

-  79  -  I 


AIR  temperature  -  °C 


FIG.  43:  EFFECT  OF  AIR  TEMPERATURE  AND  RELATIVE  WIND  SPEED  ON 
REPORTED  ICING  SEVERITY  FOR  FISHING  VESSELS  STEAMING 
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FIG.  44:  EFFECT  OF  AIR  TEMPERATURE  AND  RELATIVE  WIND  SPEED  ON 
REPORTED  ICING  SEVERITY  FOR  FISHING  VESSELS  AT  LOW  SPEEDS 

(4  KNOTS  AND  BELOW) 
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FIG.  45;  MERTINS'  CHARTS  OF  ICING  ON  FISHING  VESSELS  AT  SLOW  SPEEDS 
IN  WINDS  OF  BEAUFORT  FORCE  6  12 

(REPRODUCED  f  ROM  Rt  F  4bl 
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FIG.  47:  COMPARISON  BETWEEN  CALCULATED  RATE  OF  ICING  AND  THAT  ESTIMATED 

FROM  TRAWLER  ICING  REPORTS 
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SIMPLII  IKI)  K)KMM,AK  K)H  TOT \L  COIJ.I.C  I'lON  1,1  I  ICII  NCV 

oi  (:vlim)i:ks  and  hkctancu  i,ai{  hodii.s 


In  an  attempt  to  simplify  the  calculation  of  the  droplet  collection  efficiency  of  cyliiKler.s  ami 
rectangular  bodies  and  to  obviate  the  use  of  curves  such  as  presented  in  References  14  and  19  the  fol¬ 
lowing  formulae  are  proposed.  These  formulae  have  been  dt'rived  for  thi-  rather  limited  range  of 
atmospheric  conditions  appropriate  to  trawler  icing  and  the  icing  of  static  structures. 

The  properties  of  air  at  -  5°C  and  100  kPa  (750  mm  Ilg)  have  been  used,  and  the  deviations 
from  these  conditions  that  are  to  be  expected  in  this  context  have  negligible  effect  on  the  final  result. 
Thus  only  the  variables  of  velocity,  droplet  diameter,  and  cylinder  diameter  (or  body  width )  arc  in¬ 
volved  in  the  formulae'  and  are  combined  into  a  single  parameter,  f  (zeta),  of  the  form: 

f  =  V  d‘i  D'  (A.l) 


when*  V  =  velocity 

d  =  droplet  diameter 

1)  =  cylinder  diamet(*r  (or  body  width) 

'I'Ik'  indices  p,  q  and  r  must  be  chosen  such  that  the  families  of  curves,  resulting  from  the  use 
of  the  K  and  0  (or  Re^, )  parameters  as  used  in  Reft'rences  14  and  49  all  collapse  to  a  single  cunc  w  hen 
the  collection  efficiency,  ,  is  [)lott<'<i  against  f. 

It  was  found  that  this  n'ciuircment  was  .satisfied  over  a  representative  range  of  values  of  d 
and  I)  (with  the  exception  of  the  simultaneous  combination  of  vi'ry  small  values  of  all  three  fiarameters ) 
by  the  following  values  of  the  indices; 

p  =  ().(') 

(1  =  i.b 
r  =  1 

V"  '■  ,1'  '■ 

whence  C  l.\,la) 

I) 


I  he  units  used  throughout  m  evaluation  f  are: 


for  V’  m/s,  for  d  pm,  and  for  I)  m 


I  o  conform  ' o  the  obvious  refpiiremcnl.s  that  when  \'  or  d  is  very  large,  or  w  hen 
small  (i.e.  when  f  ‘  I  he  collect  ion  efficiency.  K,„  ,  approaches  unitv,  and  that 
f  •  0,  a  relationship  rietween  f  and  of  the  following  form  was  adopted 


I )  IS  very 
0  when 
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K„,  =  -  ^  -  for  f  ^ 

c\ 

1  +  “  “ - 

(f  P’)" 


K 


in 


0 


for  f  - ;  p 


I  A. 2  I 


wluTt'  o,  p  aiKi  n  ;ir(“  cotistants. 

Cylinder  Collection  Kffieiency 

Kvaluation  of  a,  ^  and  n  for  the  cylinder  eolh'etion  effli-iency  data  prcscntcci  in  KcfiTc  iice  1  1 
leads  to  the  follovvinfj  ri'lations: 


1  + 


1 

8200 
f  3200 


f  3200 
f  +  27000  ’ 


3200 


0 


S  •  3200 


(A, 3  I 


This  curve  is  |)res(‘nted  in  Fi>;ur(>  A.l,  to^tether  with  actual  values  of  ly,  ev  aluated  usiiu; 
Hef(>renc('  11  ovt'r  the  following'  ranue  of  values: 

\’  from  3  m/s  to  (50  m/s 
d  from  20  pm  to  1000  pm 
I )  from  0.03  m  to  1  m 

•As  the  points  indicate,  the  vakus  of  are  collapsed  w('ll  hy  IA|ual  ion  i  .  \ .  1  a  |  onto  a  si n pic 
curve,  w'itli  tile  exce[)tion  of  a  few  points  wlucii  (ie  lielow  tiie  curve,  p.irt  ii'ufarlv'  at  its  lower  end,  and 
w  hich  correspond  to  comhiiuit ions  of  the  lower  vtilues  of  ilm  parameters  \'.  d  and  I ). 

The  curve  piven  liy  lAiuation  (.\.3)  is  seen  to  correspond  to  hettcr  than  0,0  1  in  k,,,  value  to 
the  mam  proup  of  points,  while  the  e.x'rancous  fioints  he  no  more  than  O.OS  m  1'.^,  value  from  the 
cun  I'. 


Kecttinpular  Body  Collection  efficiency 


I'he  rcctanpular  hody  collection  cffiiaency  data  of  Ueference  pl  vvcri'  used,  and  led  to  the 
follovvinp  relations: 


1  S  2S00 

1  l.htIO  V'  '  I  1  TOO  ' 

1  ♦ 

V  2,soo 


v'  •  2800 


( .\,  1  1 


2800 


I'lpUl  e  ,  \  ,  I  i)re>cMI  s  t  ills  Cline  topel  her  Wit  h  v.lllles  of  K.  cv  aliiat  ed  minp  Iv  e  fere  111  e  p.  1  ov  er 
the  ^amc  ranpe  of  the  v  ana  hies  an  d  as  for  the  cvliiidcr  data,  however  the  li...  ]  v  w  idl  h  1 1 » l  w  ,e 
e.\t  ended  to  3in  .  The  ly.^  data  eol  lapses  onto  a  sinple  eiin  e  rat  her  het  t  er  t  li,m  lor  tile  ey  I  inder  ea  -  e, 
except  will  11  the  droplet  HeynoliU  .\umlier  I  lie  )  is  less  than  alioiit  20.  in  wha  li  ease  Ihe  points  fall 
liel,  lu  the  mam  proupmp  of  points,  Ihe  eiine  i  K(|,  i  ,\.  1  1 1  .iprei's  well  v\  ii  ii  the  m.nii  In  idv  of  p.  mil  - 
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APPKNDIX  l{ 

DERIVATION  OF  FXPRF.SSION  FOR  PRKDJCJ’JNO 
THK  ICINO  H  ATH  OF  FISIIINCi  TRAW  LFKS 

B.1.0  HEAT  BALANCE  AT  THE  ICING  SURFACE 

Taking  only  the  primary  heat  transfer  processes  acting  at  the  icing  surface  into  consideration, 
and  assuming  that  the  ice  formation  is  a  continuous  steady-state  process,  an  e(juilil)rium  heat  balance 
at  the  icing  surface  may  be  formulated  as  follows: 

q,  +  +  q^.  +  q,.  +  =  0  (B.li 

where  q,  =  latent  heat  of  freezing  a  certain  fraction  of  the  impinging  water. 

q^^  =  heating  (or  cooling)  of  the  impinging  water  to  the  ecjuililirium  surface 

tem[)erature, 

q^  =  heat  transfer  l)y  conv(>ction  to  the  surrounding  air, 

q^,  =  heat  transfer  by  eva))oration  to  the  surrounding  air. 

q  I  =  heal  conducted  to  or  from  the  icing  surfat’c  through  lh(>  underlying 

structure. 

Each  of  these  terms  will  initially  be  considered  individually  in  the  following  suh-sections, 
B.1.1  I'he  Latent  Heat  I'erm,  q, 

The  rate  of  water  im[)ingement  on  unit  area  of  a  surface  may  la*  ('x))ressed  as: 

=  EwV  (H.2) 

where  K  is  the  ciilh'ction  efficicmcy, 

w  Is  the  litjuid  water  content  of  the  air, 

V  IS  the  velocity  of  the  air  relative  to  tile  object  or  surface. 

I'he  CO  I  lei' 1 11)11  efficiency  e.xfiresse.s  t  he  degree  to  which  the  water  dio)  is  m  the  air  ini)  >muc  on 
an  object  in  their  |iath.  anil  are  not  deflected  hv  the  airllow  around  the  object .  The  eolleel  ion  et  f’eie'ii  \ 
in(‘rease>  with  rclat  iie  velocil  \'.  and  with  dro|)  si, a',  and  ilecrcases  with  the  si/.e  ol'  !  he  ohier  i  on  w  li  n  'h 
the  (lro[  o  impinge.  Hi ‘cause  the  size  of  t  hi'  w  ater  drops  involved  in  sle  i  lemu  due  i  o  sea  spray  is  hn  gi 
I  1  nun  ).  they  will  he  denect  ed  little  and  the  collect  ion  efficiency  will  lie  assn  n  led  to  In  1  III) 

The  I  lisi  III  ml  loll  of  I  Irop  sizes  in  the  s)iray  is  unknow  n  and  vanes  with  wind  -pi  ii  j ,  -i ;  p  pi  eo 

t  V  I  le  and  si/e  of  \  e,,sel .  !  he  local  loll  on  the  vessel,  and  with  height  above  the  sea  siirf.ii  I  I  )a  I  a  o yn  gei  I 

I  iv  Borisenkov  (  l\  ef.  .'I  I  suggests  I  he  I  Irop  diamei  i-r  vanes  fri  mi  1  .0  to  d..'i  iiini  1 1  In  ,i\  eiage  i-  p.  1  nun  i  i.. 
a  re-n  h  o  f  spray  inu  fri  mi  -well  w  .ives  I  low  ever,  in  the  I  mw  of  t  lie  sliii '  the  w  .iter  will  In  eni  onni  i  i,  -  j  in 
I  'l  .e  form  of  ■  ■'gllsllis''  or  "ra'iM.  ies"  for  W  llleli  die  eoni  epi  I  if  drop  s|..-,'  -eeln-  lllappl'i  iprial  '  el  !  1 1,  ■!,  . 

lor  dll'  wall!  of  .1  i  let!  er  coni  e(  a  .  It  Will  lie  eniiiloyeil  m  |  pi  i-  anal  v  -i  -. 

I  he  I  ii  |U  id  'w  at  er  I  oiii  I 'lU  1  if  I  he  ray  clo.id  I- eiju.ill  V  ihffiriili  t  o  i  pi.int  1 1  >  in.inv  'i;na|ni  w  ,i'. 

how  ev  I  r.  Kaeliiinn  I  B  ef  1  7  I  -iiggest  s  t  hat  waveheiglii  i-  i  In- ehief  hn  i  or  gov  erninr  :  h  i  p.namitii  \ 

w  It  h  drop  -I  /  e.  It  ni.iv  lie  a,--'  miei  I  l  li,il  die  vv  .it  i  r  rout  nii  o  al.-o  fin  let  ion  of  da  si./i  .  I  m  -  and  mi  d 
of  dll'  w‘--e|,  and  .iPo  .1  f  mil  I  ion  ol  heigh i  .li mve  the  -e.i  siirf.ii  e.  K.n  bin  in  propo-e  ;  P.u  | ,  ,r  -.ri i m 
-peei  1  .  o  t  li  s  1.  n ,  i| ,.  and  I  a  '.r  I  mg-  o  I  I )  '  in  .  a  |  n'oi  mri  aaial  rel.ii  lon-i  np  oe  ,is.i  n  m 'd 


r  .vivCiguieiVi  BLAiNE -ivOT  ftjLj'iLL 
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w  =  10’’  kg/m'  (Ik.'ii 

where  is  the  wave  height  in  metres. 

However,  in  the  icing  model  developed  here,  it  was  found  that  a  liquid  water  content  one 
sixth  of  that  given  by  expression  (B.3)  resulted  in  a  significant  improvement  in  correlation  witli  observed 
results.  .As  a  result,  the  expression  for  liquid  water  content 

w  =  1.7  X  lO-*  kg/m'  (B.3a) 


has  been  adopted. 


If  a  fraction,  n,  of  the  impinging  water  freezes  on  the  surface,  then  the  rate  of  ice  growth  |ier 
unit  area  i.s: 


Ri  =  nR^^  =  nKwV 

and  if  t,  is  tlu'  latent  heat  of  freezing,  then  the  heat  released  is: 

9,  = 


IH.  I  I 


iB.dl 


I'he  .Sensible  Heat  Term, 

If  the  equilllirittm  temperature,  t.^.  of  the  icing  surface  on  which  the  water  drops  impinge  is 
assumed  to  be  the  etiuilibrium  freezing  temperature,  t,.of  water  of  the  appropriate  salinity  li.e.  ignoring 
any  temperature  gradient  ai  ross  any  water  film  on  the  icing  surface),  aiul  tlu'  drop  temperature  imme¬ 
diately  prior  to  impingement  is  t  j .  thim  the  heat  given  to  the  imi>inging  water  is: 

‘lu  ~  R,v  '*a  -  Oktii 

wiiere  Is  tlv'  specific  heat  of  the  impinging  water  >  101)0.]  kgC 

for  s  ■-  (•'  It)  C. 

It  should  not  lie  overlooked  liowever  that  some  of  the  impinging  water  w  ill  likely  splash  or 
run  off  the  icing  surface  before  reaching  the  e.uiililirium  temperature,  t^.  Thus  a  splash  factor,  s.  sho.ild 
ap(irofiriatel>'  be  introduced  mt o  l'l((ual  ion  (  B.b )  as  follows; 

((,,  (1  s)  R,^  c,,  (1,1  tj  I  li.ba  1 

I  low  ever,  there  seems  lit  I  le  way  of  knowing  (he  magnil  ude  of  the  s,  but  .li.ielinrn' 

I  Ref .  17  I  suggest s  t hat  about  0.7.7  of  the  water  film  runs  off  without  1  reez.mg  and  assumes  there  is  no 
contribution  whatever  to  the  heal  lialanee  from  this  water.  In  f.ict.  it  is  likely  lliat  the  lemperature  of 
the  w  ater  film  is  close  to  t  liat  of  the  lee  surface,  and  t  hat  it  is  only  the  w  a' er  that  splashes  and  is  blow  n 
away  at  the  moment  of  im|nngement  that  takes  no  part  in  tne  thermal  process. 


Becaioe  of  these  uncertainties  the  form  -if  Kqiialion  (  B.bl  will  be  enqiliw  ed  n,  the  follow  im, 

analysis. 


I  he  tenq  lerat  lire  of  the  inioniging  water  dro|)s,  t  , ,  i-  a  In  mi  ion  of  man  \  \  .n  i,  hies  indui  I  ii  i; 
t  hi-  ^ea  surface  t  emperat  nre,  the  air  lemperal  nre  I  he  drop  size  and  I  he  t  inie  ot  fl  igh  l . 
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An  expression  for  the  drop  temperature  may  be  derived  by  following  the  procedure  due  to 
Hardy  (Ref.  50). 


The  rate  of  change  of  temperature  of  the  drop  due  to  convection  is: 


Nu 


(ta-  t) 


where  t 

T 

Nu 

K 

Pw 

Cw 

d 


instantaneous  drop  temperature 

air  temperature 

time 

Nusselt  No. 

thermal  conductivity  of  air 
density  of  water 
specific  heat  of  water 
drop  diameter 


(B.7) 


The  rate  of  evaporation  is  given  by: 


where  m 


D 

Pa 

Mw/Ma 
e.  Ca 
p 


dm  K-  e\ 

—  -  Na.dD,.  j  ( 

mass 

diffusivity  of  water  vapour  in  air 
density  of  air 

ratio  of  molecular  weights  of  water  and  air  (0.622) 

vapour  pressure  at  surface  of  drop,  and  in  surrounding  air  respectively 

barometric  pressure 


But  the  temperature  change,  dt,  resulting  from  an  evaporative  mass  loss,  dm,  from  a  drop  of 
diameter  d  is: 


dt 


6  fi. 


Pw 


dm 


(B.9) 


where  Cy  =  latent  heat  of  vaporization  (2.5  X  10*  J/kg). 
Thus,  the  rate  of  change  of  temperature  due  to  evaporation  is: 


egyPPa  /Ca  -  e\ 

M„  \  p  y 


(B.IO) 


Combining  Expressions  (B.7)  and  (B.IO)  to  obtain  the  overall  rate  of  temperature  change: 
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dt  _  6Nuk, 

dr 

6  Nu  kj 

"  p  c  ■  ■" 

where  X,  is  a  factor  to  allow  for  evaporation: 

®v  ®a  ~  ® 

X,  =  1  +  0.622  -  - 

PCp  ta  -  t 


(ta-  t)U  + 


M^C^Dpj 

MjkaP 


f;^)l 


V  /J.  iV 


(B.ll) 


(B.lla) 


(B.12) 


k 

since  D  - 

PaCp 


and  where  Cp  =  specific  heat  of  air  (1.005  X  10^  J/kgK). 


Since  X,  is  a  weak  function  of  temperature.  Expression  (B.lla)  may  be  integrated  assuming 
X,  is  a  constant,  resulting  in: 


(t^  -  td )  6  Nu  ka 


(B.13) 


or  td  =  t^  +  (t„  -  tj  e 


/  6  Nu  \ 

\  ‘  / 


(B.14) 


where  td  is  the  water  drop  temperature  at  time  r. 


To  test  the  validity  of  the  assumption  that  the  factor  X,  may  be  treated  as  a  constant  in  the 
integration,  results  obtained  using  Equation  (B.14)  are  compared  in  Figure  B.l  with  the  results  of  a 
numerical  integration  of  Equation  (B.ll).  Also  shown  for  comparison  is  the  effect  of  ignoring  the 
evaporative  effect  (convection  only).  It  is  clear  that  the  assumption  of  a  constant  X,  does  not  result 
in  any  serious  errors  in  the  computation  of  drop  temperature. 

The  appropriate  time  or  distance  of  flight  of  a  drop  before  striking  some  part  of  the  ship  is 
uncertain,  as  is  also  a  representative  drop  diameter.  Clearly,  larger  drops  will  stay  airborne  for  shorter 
times  than  smaller  drops,  so  that,  together  with  the  increase  in  the  value  of  the  time  constant 

I X,^  in  Equation  (B.14)  resulting  from  a  larger  d,  larger  drops  experience  consideralily  le'ss 

cooling  than  smaller  drops,  a  circumstance  which  results  at  times  in  the  decks  and  lower  parts  of  a  vessel 
experiencing  little  ice  build-up  while  appreciable  ice  forms  on  rigging  and  superstructure. 


,6  k,  Nu 


For  the  sake  of  this  analysis,  however,  an  effective  drop  diameter  of  2  mm  is  assumed.  If  is 
further  assumed  that  the  effective  relative  velocity  between  the  drop  and  its  surrounding  air  is  the  fall 
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velocity  of  the  drop,  which  for  a  2  mm  drop  is  6.5  m/s  (see  Ref.  51 ).  Thus,  using  the  expression  for 
the  convective  heat  transfer  for  a  sphere  given  by  Mc.\dams  (Ref.  52)  i.e.  Nu  =  0.37  Re®  *,  the  Nusselt 
number  is  23  and  hence  the  time  constant  is  approximately  5/X,  seconds. 

The  flight  time  of  the  drops,  t,  was  expressed  in  Reference  42  as  x/V  (i.e.  the  distance  carried 
by  wind  of  velocity  V)  and  it  was  suggested  that  x  might  be  related  to  vessel  length  as  a  means  of  intro¬ 
ducing  vessel  size  into  the  analysis  however,  this  is  perhaps  not  the  appropriate  place  to  introduce  size 
since  the  effect  would  be  of  increasing  icing  rate  with  increasing  vessel  length,  a  situation  that  does  not 
always  accord  with  observation.  Perhaps  some  modification  to  the  liquid  water  content  expression 
(Eq.  (B.3))  would  be  the  more  appropriate  way  of  introducing  both  the  size  and  lines  of  the  vessel. 

More  knowledge  of  the  effect  of  size  and  lines  on  the  rate  of  icing  is  first  needed. 

In  this  analysis  the  flight  distance,  x,  is  assumed  constant  at  20  m,  whence  the  effective  flight 
time  of  the  drops,  t,  becomes  20/V,  a  figure  that  results  in  the  best  correlation  between  the  analytical 
and  observed  rates  of  icing. 

The  surface  temperature,  tj,  was  stated  earlier  to  be  the  equilibrium  freezing  temperature  of 
water  of  the  appropriate  salinity.  As  ice  forms  at  the  surface,  salt  is  rejected  from  the  ice  crystal  lattice 
so  formed,  and  the  water  film  on  the  surface  of  the  ice  manifests  an  enriched  salt  content.  Since  a 
definite  relationship  exists  between  the  salinity  and  the  freezing  temperature  of  the  brine,  it  follows 
that  the  surface  temperature,  t,,  must  be  that  appropriate  to  the  salinity  of  the  brine  on  the  surface  of 
the  ice,  and  not  that  appropriate  to  the  parent  sea  water,  tf.  Further,  Tabata  (Ref.  2)  has  shown  that  the 
chlorinity  and  freezing  temperature  of  the  brine  run-off  is  a  function  of  the  freezing  fraction,  n(=  Rj/R^  ). 
In  fact  his  data  suggests  the  simple  relationship; 

t5  =  (l+n)tf  (B.15) 

although  this  likely  does  not  hold  as  n  1 . 

The  freezing  temperature  is  related  to  the  salinity  by : 

tf  =  -0.002  -  0.0524S-  6.00  X  lO'^S^  (B.16) 

where  S  is  in  parts  per  thousand. 

This  expression  is  within  0.002‘’C  of  the  values  given  in  Reference  39  for  salinities  between  0%o  and  40%o . 

Thus,  incorporating  Equations  (B.2)  and  (B.3a)  in  Equation  (B.6): 

q,^,  =  1.7  X  10-’  H„c,„  V  (tj  -  tj)  (B.6b) 

where  t^  and  t,  are  given  by  Equations  (B.14),  (B.15)  and  (B.16). 

B.1.3  The  Convective  Heat  Term,  q^ 

The  convective  heat  transfer  is  expressed  in  terms  of  a  heat  transfer  coefficient,  h,  and  the 
temperature  difference  between  the  surface  and  the  surrounding  air; 

q,  =  h  (t3  -  t,)  (B.17) 


In  determining  an  expression  for  the  convective  heat  transfer  coefficient,  a  certain  amount  of 
approximation  is  necessary  since  the  shape  and  size  of  the  icing  surface  is  undefined. 
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Por  a  flat  plate  in  turbulent  flow  parallel  to  its  surface  (c.f.  the  deck  of  a  vessel)  an  average 
heat  transfer  coefficient  over  a  length  L  is  given  (Ref.  53)  by; 

h  =  0.037  —  Pr‘/3Re®® 

L 

If  a  mean  “film”  or  boundary  layer  air  temperature  of 
reduces  to: 

VO.8 

h  =  6.44  -  W/m^K  (B.19) 

LO.2 

In  the  case  of  a  cylinder  of  diameter  D  in  crossflow  at  a  Reynolds  number  of  between  50,000 
and  200,000  the  heat  transfer  coefficient  may  be  expressed  (Ref.  54)  by: 

k^ 

h  =  0.0208  —  Re"*  I**  (B.20) 

D 

An  expression  that  provides  the  same  index  to  the  Reynolds  No.  as  in  Equation  (B.18)  and  is 
within  1%  of  Equation  (B.20)  over  its  applicable  Reynolds  number  range  is: 

^a 

h  =  0.0244  —  Re"  *  (B.21) 

which  at  a  mean  film  temperature  of  -  S^C  reduces  to; 

VO.8 

h  =  4.75  -  W/m^K  (B.22) 

D".2 

Expressions  (B.20)  through  (B.22)  are  appropriate  to  smooth  cylinders  in  air  flows  of  low 
turbulence.  The  effect  of  roughness  and  of  free-stream  turbulence  can  be  considerable  (Ref.  54),  and 
since  the  air  flow  over  a  ship  is  highly  turbulent,  it  does  not  seem  unreasonable  to  increase  the  coeffi¬ 
cient  of  Equation  (B.22)  arbitrarily  and  equate  it  to  that  of  Equation  (B.19).  If  then  a  characteristic 
dimension  is  arbitrarily  chosen  as  3  m  in  each  case,  there  results  the  expression  for  both  flat  and  cylin¬ 
drical  surfaces; 

h  =  5.17  V"*  V/m^K  (B.23) 

Because  of  the  power  (0.2)  to  which  the  characteristic  dimension  is  raised,  the  heat  transfer 
coefficient,  h,  is  not  very  sensitive  to  this  dimension,  being  within  10%  of  its  nominal  value  given  by 
Equation  (B.23)  for  values  of  L  from  1.8  m  to  5  m. 

Equation  (B.17)  can  finally  be  written; 

q,  =  5.17  V"*  (t,  m  (B.24) 


(B.18) 

-  5°C  is  assumed,  this  expression 


where  t,  is  the  air  temperature,  and 

t,  is  the  surface  temperature  defined  by  Equation  ( B.  1 5 ). 
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B.1.4  The  Evaporative  Heat  Term,  q. 


The  evaporative  heat  transfer,  ,  can  be  evaluated  by  analogy  with  the  convective  heat 
transfer  (Ref.  52): 


9e 


e  Cy 
PCp 


-  es) 


(B.25) 


where  Pr  = 
Sc  = 
e  = 


P  = 
Cv  = 


Prandtl  number  =  0.711 
Schmidt  number  =  0.595 

ratio  of  molecular  weights  of  water  vapour  and  dry  air  (0.622) 
atmospheric  air  pressure 

latent  heat  of  vaporization  of  water  (2.5  X  10*  J/kg) 
specific  heat  of  dry  air  (1.005  X  10*  J/kgK) 

saturation  vapour  pressure  of  moist  air  at  temperature  t^,  t^  respectively. 


Evaluation  of 


,  e,  Cy  and  Cp  results  in: 


q,  =  1731  ^  (e,  -  e,)  (B.26) 

and  if,  as  before,  h  =  5.17  V®  ®W/m*K 
and  assuming  p  =  100  kPa 
then: 

qe  =  89.5  (63  -  Cj)  V"  *  W/m*  (B.27) 

where  ej  and  ej  are  in  kPa,  and  V  is  in  m/s. 

The  saturation  vapour  pressure,  e,  may  be  expressed  in  kPa  as  a  function  of  temperature,  t(°C), 
by: 

e  =  ^((1.9226  X  lO-’t  +  2.4545  X  10'’ )t  +  1.4224  X  lO'*)  t  +  0.044436^1  +  0.61094  (B.28) 

This  4th  order  polynomial  regression  equation  is  within  '/4%  of  the  values  given  in  Reference 
39  in  the  temperature  range  -25°C  to  +10°C. 


B.1.5  The  Conductive  Heat  Term,  q. 


In  this  analysis  the  heat  conducted,  qg,  to  or  from  the  icing  surface  through  the  ice  layer  and 
the  underlying  structure  will  be  assumed  negligible,  although  frequently  instances  may  be  observed  on 
iced  ships  of  inadequately  insulated  crew  spaces,  as  indicated  by  a  lesser  quantity  of  ice  in  these  areas, 
thus  clearly  indicating  that  in  these  instances  heat  is  being  conducted  through  the  structure. 


-96- 


B.2.0  The  Total  Heat  Balance  and  Rate  of  Icing 

The  overall  beat  balance  (Eq.  (B.l))  may  now  be  written  using  Equations  (B.5),  (B.6b), 

(B.24)  and  (B.27); 

CfRj  +  0.68H^V(td  -tj  +  5.17V°8  j(t^  -  +  17.3  (e^  -  e^)}  =  0  (B.29) 

and  evaluating  C,  as  3.33  X  10®  J/kg,  the  rate  of  icing  is  given  by: 

Rj  =  2.04  X  10-«  H^V(ts  -  td)  +  1.55  X  10*®  V"  *  -  t^)  +  17.3  (e^  -  ejjkg/m^s  (B.30) 

Or,  alternatively  in  terms  of  rate  of  growth  of  ice  thickness,  Nj.- 

Rj 

Nj  =  —  =  2.3  X  10-®  V(t,  -  tj)  +  1.74  X  lO'*  V"-®  ((t^  -  t J  +  17.3  (e,  -  e3))m/s  (B.31a) 
where  the  density  of  ice,  pj,  is  assumed  to  be  890  kg/m* . 

If  the  rate  is  to  be  expressed  in  the  more  convenient  units  of  mm/h,  then: 

Nj*  »  8.26  X  10-*  H^V(t,  -  td)  +  6.28  X  10'*  V"  ®  {(t^  -  tj  +  17.3  (e^  -  ej)mm/h  (B.31b) 

In  these  equations  the  wave  height,  ,  is  in  metres,  the  relative  wind  velocity,  V,  is  in  m/s, 
the  tempe^atures  are  in  °C  and  the  vapour  pressures  are  in  kPa. 

If,  as  is  more  usual,  the  velocity  is  given  in  knots.  Equation  (B.31b)  becomes: 

Nj*  =  4.25  X  10-*  Vk,(t5  -  td)  +  3.69  X  10'*  Vk,®®  {(t^  -  t^)  +  17.3  (e^  -  e3)|mm/h  (B.31c) 

To  test  this  model  of  icing  under  conditions  of  sea  spray,  the  icing  rate  as  calculated  by  Equa¬ 
tion  (B.31c)  using  the  conditions  applicable  to  39  reported  icing  incidents  (Table  VIII)  has  been  com¬ 
pared  in  Figure  43  with  the  mean  icing  rates  derived  from  these  reports. 


B.3.0  Practical  Implementation  of  the  Rate  of  Icing  Equation 

In  evaluating  those  expressions  that  depend  on  the  vapour  pressure  of  the  atmosphere  (i.e. 
Eq.  (B.12),  the  evaporation  factor  in  the  drop  cooling  equation,  and  Eq.  (B.27),  the  evaporative  heat 
term  in  the  heat  balance  at  the  icing  surface),  it  has  been  assumed  that  the  air  is  90%  saturated  rather 
than  saturated,  so  that  in  both  Equations  (B.12)  and  (B.27)  the  term  is  replaced  by  the  expression 
0.9  e,. 


Because,  from  Equation  (B.15),  the  effective  surface  temperature,  t^,  is  a  function  of  the 
freezing  fraction,  n(=Ri/R^^, )  and  hence  a  function  of  the,  as  yet  unknown,  icing  rate,  it  is  necessary 
to  employ  a  simple  iterative  procedure  to  determine  the  appropriate  value  of  t^  to  use  in  Equations  (B.6), 
(B.24)  and  (B.28).  The  initial  value  of  the  freezing  fraction,  n,  was  assumed  to  be  zero,  i.e.  t^  =  t,  .  Sub¬ 
sequent  iterations  used  the  computed  vjdue  of  Rj  (Eq.  (B.30))  to  determine  a  new  value  of  n,  until  suc¬ 
cessive  values  of  n  differed  by  less  than  an  arbitrarj'  amount  (0.0001  was  used).  If  n  assumed  a  negative 
value,  it  was  set  at  zero  and  no  icing  was  inferred. 
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Similarly,  in  determining  the  final  drop  temperature,  tj ,  it  was  initially  set  at  the  sea-water 
temperature,  t^„ ,  and  a  value  of  X^  (Eq.  (B.12))  computed;  this  was  in  turm  employed  in  Equation  (B.14) 
to  derive  a  new  value  of  tj  to  be  used  in  Equation  (B.12).  This  procedure  was  repeated  until  successive 
values  of  tj  differed  by  less  than  an  arbitrary  amount  (again  0.0001°C  was  used). 


I 

J 


r 


DROP  TEMPERATURE 
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FIG.  B.1;  COOLING  OF  2  MM  DIAMETER  WATER  DROPS  FALLING  AT  TERMINAL 
VELOCITY  THROUGH  AIR  OF  TEMPERATURE  -  10°C 
INITIAL  DROP  TEMPERATURE  +2°C 
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